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1.0 INTRODUCTION AND SUMMARY 


TRW Defense and Space Systems has performed a series of Application 
Studies for NASA Goddard Space Flight Center under contract NAS5-21188. 
These studies have focused upon three-axis attitude control system 
designs employing momentum/reaction wheels with magnetic or thruster momen- 
tum unloading. The studies encompass conceptual design, system tradeoff, 
control law development, performance analysis and simulation, and hardware 
implementation considerations. The current study is directed toward 
design and performance evaluation of the LANDSAT D attitude control sys- 
tem (ACS). Control and configuration of the gimballed Ku-band antenna 
system for communication with the Tracking and Data Relay Satellite (TDRS), 
and control of the solar array drive are considered a part of the ACS and 
are also addressed here. 

This report is an Interim Technical Report and documents only some 
aspects of the study. A complete final report will be published at the 
end of the study and the present material will then also be included. 

The LANDSAT D mission, a follow-on to LANDSAT's A, B, and C, is cur- 
rently in the planning stages and it will utilize the NASA Multi-Mission 
Modular Spacecraft. It will fly at an altitude of 705 km in a sun syn- 
chronous, 9:30 a.m. (descending node) orbit. Complete earth coverage is 
obtained in 16 days; exact revisit time is 233 orbits. Launch is tenta- 
tively planned for early in 1981. The payload consists of a Thematic 
Mapper (TM) and a Multi -Spectral Scanner (MSS). Both sensors collect 
spectral radiation data of (lighted) land masses. This data will be used 
to produce land use theme maps to better manage the earth resources. 
Potential users have a broad spectrum of interests including agriculture, 
mineralogy, and urban planning, to cite a few. The spatial resolution of 
the Thematic Mapper is 30 meters corresponding to the 30 meter square 
footprint on the ground of a square detector in the instrument. The 
resolution of the MSS is 80 meters. The spacecraft pointing accuracy is 
specified as 0.01 deg (lo) witha jitter stability of +10.5pRad (lo) over 
20 minutes. The ultimately desired temporal registration (essentially 
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repeatability of ground track after 16 days) of the Thematic Mapper images 
is to be within 3 meters (la). 

The primary objective of this study is to assess the adequacy of the 
Multi-Mission Modular Spacecraft (the MMS attitude control system in par- 
ticular) to perform the LANDSAT D mission, and to identify and configure 
mission peculiar equipment that may be required. Examples are the gim- 
balled Ku-band antenna and the single sided solar array and drive. To 
this end, preliminary design studies and performance assessments were con- 
ducted using the MMS as the baseline system concept. 

Section 2 describes the LANDSAT D mission in more detail. Important 
currently available mission data is sutrenarized and a brief description of 
the MMS and the two payload sensors is provided. 

The gimballed Ku-band antenna system for communication with TDRS is 
discussed in Section 3. By means of an error analysis it is demonstrated 
that the antenna cannot be open loop pointed to TDRS by an on-board pro- 
grammer but that an autotrack system is required. After some tradeoffs, 
a two-axis, azimuth-elevation type gimbal configuration is recommended for 
the antenna. It is shown that gimbal lock only occurs when LANDSAT D is 
over water where a temporary loss of the communication link to TDRS is of 
no consequence. A preliminary gimbal control system design is also pre- 
sented. 

-A digital computer program was written that computes antenna gimbal 
angle profiles, assesses percent antenna beam interference with the solar 
array (as a function of antenna mast length), and determines whether the 
spacecraft is over land or water, a lighted earth or a dark earth, and 
whether the spacecraft is in eclipse. Besides checking for optical block- 
age of the antenna beam by the solar array paddle, another important appli- 
cation of this program is to determine when the spacecraft is over lighted 
land, since this is directly proportional to the TM electric power duty 
cycle. For one typical 15 orbit (24 hrs) case, the spacecraft spent only 
16.5% of the time over lighted land. 
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A solar array control strategy was developed that eliminates inter- 
ference with the Ku-band antenna beam. The array is driven at orbit rate 
as long as the spacecraft is within +45° of the ecliptic plane. When the 
spacecraft is outside this region the array is kept fixed at the orienta- 
tion it had at away from the ecliptic and this orientation is main- 
tained until the earth becomes dark. At that point the array is slewed 
again to its optimal position normal to the sun line and it continues to 
track the sun until the spacecraft goes into eclipse. This strategy 
requires the Ku-band antenna mast to be only 2.11 m (6.9 ft) in length in 
order to eliminate all interference. The key is the 37.5° cant in the 
array which permits the antenna to look past the array when it is within 
+45° of the ecliptic. Outside this region the array would rapidly 
obstruct the field-of-view of the antenna if it were not stopped as 
described above. The total power loss associated with this solar array 
control strategy is only 4% and it can be easily compensated for by 
increasing the array area by 4%. The 2.11 meter long antenna mast has 
been computed for a 4% increased array area. 

Section 4 sizes external disturbance torques acting on LANDSAT D. 
Aerodynamic, gravity gradient, magnetic, and solar pressur e di sturbance 
torques are considered. Ku-band antenna mast length and antenna orienta- 
tion were observed to affect the magnitude of the disturbance torques only 
insignificantly. 

The disturbance torque sizing results were then used to assess 
whether the four 20 Nms (15 ft-lb-sec) reaction wheels and the three ortho- 
gonal 100,000 pole-cm magnetic torque bars provided in the MMS ACS module 
have sufficient capability for the LANDSAT application. The conclusions 
were that the spacecraft can operate effectively in its 705 km orbit 
environment even when using only the 50,000 pole-cm windings on the mag- 
netic bars. Momentum excursions in the wheels were observed within 
+2.0|Nms of nominal and since the reaction wheel momentum storage capacity 
is 20 Nms, it is recorrmended that all four wheels be run at a biased 
speed. This avoids wheel speed reversal and associated attitude transients 
due to nonlinear wheel friction. The four wheels are speed biased such 
that a net zero momentum system results (not possible with only three 
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wheels). Should one wheel fails then the three remaining wheels must be 
operated about zero speed and the wheel speed reversal attitude transients 
of about 90 yRad (18.5 arc-sec) must be accepted. This would occur twice 
per orbit in roll and twice per orbit in yaw. 

Initialization of the fine control mode is discussed in Section 5. 

Two methods to perform fine stellar acquisition, i.e., to initialize the 
on-board stellar-inertial attitude reference system,are presented. Before 
entering the stellar acquisition mode the spacecraft is sun pointing (in 
the fine sun acquisition mode) and only a relatively coarse determination 
of the spacecraft orientation about the sun line is needed for correct star 
identification. One method uses magnetic field measurements and is cap- 
able of providing attitude knowledge about the sun line within about 2°. 

The other method correlates observed stars in a swath about the sun l,ine with 
star catalog entries and determines attitude. about the sun line to within 
approximately 0.2 degree (uncalibrated gyros assumed). The magnetic field 
measurement method is simpler and the 2° attitude knowledge is sufficient 
for unambiguously identifying the first few stars when the stellar-inertial 
attitude reference system is initialized. Using the real star field it 
is demonstrated that the linearized (extended) Kalman filt er al gorithm 
assumed to be used by the on-board attitude reference system, is capable 
of converging rapidly from large (2°) initial attitude uncertainties. 

The converged attitude reference filter predicted accuracies in the 5 to 
10 arc-sec (la) range and the gyro biases were calibrated to within about 
0.0013 deg/hr. 

The final report will address normal on-orbit ACS design and per- 
formance. The effects of flexible and moving appendages such as the Ku-band 

antenna and the solar array are also considered. Stationkeeping, autono- 
mous Thematic Mapper boresight registration and on-board software sizing 
will also be addressed. 
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2.0 LANDSAT-D -MISSION 


LANDSAT-D is the planned fourth mission in the LANDSAT program whose 
primary objective is to evaluate and monitor the earth's resources from 
satellite observations. Formerly, this program was known as the Earth Re- 
sources Technology Satellite (ERTS) program. The spacecraft to be flown 
on the LANDSAT-D mission, as well as the mission itself, have in the past 
also been referred to as LANDSAT Follow-On (LF/0). The primary mission 
payload sensor for LANDSAT-D (L-D) is the Thematic Mapper (TM), a large 
multi -spectral instrument which uses solid state detectors to produce 
earth imaging data which is digitized and transmitted to the ground in a 
120 megabit per second serial data stream. Potential users of TM data 
have a broad spectrum of interests including agriculture, mineralogy, and 
urban planning, to cite a few. 

Agriculture as it relates 'to the economy and to the world food supply 
is one of the driving forces behind the LANDSAT program. Data from the 
Thematic Mapper can be used to determine how much wheat, corn, barley, etc 
have been planted, how well the crops are doing, and combined with water- 
shed data also collected by the Thematic Mapper, predictions of yield can 
be made. The system can collect this data on a world-wide basis. The 
public benefit is expected to be enormous. Shortages can be predicted and 
additional acerage planted - perhaps in a different hemisphere. Surpluses 
in one crop may be avoided by planting different crops. Effects of dis- 
ease or insects can be assessed and hopefully minimized by early detection 

Two LANDSAT spacecraft are currently in orbit, LANDSAT's A and B, 
with an instrument known as the Multi -Spectral Scanner (MSS) being their 
primary payload sensor; the MSS is basically a smaller and less advanced/ 
less powerful version of a Thematic Mapper type instrument. LANDSAT's A 
and B are producing image data which is primarily being used for research 
in remote sensing technology, data interpretation techniques, etc. A 
third spacecraft, LANDSAT-C, is being prepared for launch with a modestly 
upgraded Multi -Spectral Scanner instrument which will include along wave 
infra red (10.5 to 12 microns) sensing capability in addition to the four 
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spectral bands carried on LANDSAT's A and B. This upgraded version of the 
MSS will also be flown as a secondary payload on LANDSAT-D. It will pro- 
vide continuing service to users who have not yet upgraded their ground 
stations and data processing capabilities to handle Thematic Mapper data. 

It will also provide an opportunity to directly compare spectral signature 
data from the two instruments, maximizing learning transfer and enhancing 
capability to compare Thematic Mapper data with historic MSS data. 

2.1 LANDSAT-D Orbital Data 

The LANDSAT-D (L-D) spacecraft will fly in a circular, 9:30 a.m., de- 
scending node, sun synchronous orbit. The terminology "sun-synchronous" 
refers to the property that the line of nodes of the orbit advances at the 
same rate as the earth’s orbital rate about the sun. The orbit altitude 
is 705 km (380 nautical miles) and the orbit plane inclination is 98.2 de- 
grees. The orbital period can be calculated as 1.648 hours (98.88 minutes), 
or stating it differently, there are exactly 14 9/16 orbits in a 24 hour 
period. Figure 2-1 shows the L-D orbit relative to the earth equatorial 
plane at vernal equinox, illustrating the 9:30 a.m. sun synchronous orbit 
feature. It simply means that throughout the year the spacecraft will 
cross the equator on its north-to-south descending pass at 9:30 a.m. local 
ground time. Since the orbit is almost polar, the ground time is close to 
9:30 a.m. over the entire sun lit side of the earth, and close to 9:30 p.m. 
over the entire dark side of the earth. 

Nominally the Thematic Mapper (TM) data is transmitted in real time 
from the L-D to a geostationary Tracking and Data Relay Satellite (TDRS) 
which retransmits the data to the ground. By 1981 when L-D is scheduled 
to be launched, there will at least be two operational TDRS's in geosyn- 
chronous equatorial orbit. They will be positioned above the equator at 
the following longitudes. 

• TDRS East: 41“ West Longitude 

• TDRS West: 171° West Longitude 

The TDRS equatorial orbit is also indicated in Figure 2-1. 
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LAND$AT-D Orbit Summary 

• 9:30 a.m. , descending node, sun synchronous orbit 

• Inclination i = 98,2 deg 

• Altitude h = 705 km 

t Orbit Period T = 1,648 hrs 

0 

• Orbit Rate = 0.06068 deg/sec 

North 



Figure 2-1. LANDSAT-D Orbit 









Solar time on earth repeats at integer multiples of 24 hours. Let 
be the L-D orbit period in hours. Then the smallest number of orbits 
n such that n is an integer multiple of 24 hours is the L-D revisit 
time, i.e., the satellite will be exactly above the same point on earth 
after exactly n orbits (ideal orbits assumed). If the north-south ground 
path illuminated by the payload sensor field of view is chosen to be of 
appropriate width, then the revisit time is also the time it, takes to ob- 
tain total earth coverage, i.e., obtain data from the entire earth surface, 
with the exception of two small regions about each of the earth's poles 
which cannot be seen because the orbit inclination is 98.2 degrees, that 
is, more than 90 degrees. Since the orbit period of L-D is 


T 


0 


24 

14 9/16 


(24) x (16) 
233 


hrs 


one must find the smallest integer n such that 


(24) X (16) 
233 


n 


24 m 


where m must be an Integer. This leads to m = 16 and n = 233. Thus the 
revisit time is 16 days during which the spacecraft completes 233 orbits. 

L-D is scheduled to be launched in 1981 on a Thor Delta 3910 booster, 
but the spacecraft and payload designs must also be compatible with future 
launch, resupply, and retrieval by the Space Shuttle. L-D mission life 
has been preliminarily specified as 3 years. 

After orbit insertion of the spacecraft, initial orbit corrections 
are performed to assure orbit- circularity and correct altitude and orbit 
inclination. Thereafter, aV stationkeeping maneuvers will only be per- 
formed to make up for orbit decay due to atmospheric drag. A + 6.5 km 
E-W ground track drift due to spacecraft altitude and attendant orbit pe- 
riod change is the maximum permissible error. Using a conservative drag 
model, this implies that the orbital velocity must be corrected by apply- 
ing a AV = 0.35 m/second at most every 16 days. 
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2.2 Spacecraft Configuration 


Because the TM Is a larger payload sensor than the MSS, the LANDSAT-D 
mission will require a new spacecraft design to support the TM's greater 
weight, provide more power, provide more accurate pointing control, and to 
relay higher data rates to the ground. It has been decided by NASA that 
the L-D mission will be performed using the Multi -Mission Modular Space- 
craft (MMS). The MMS is a concept that has been developed by NASA Goddard 
Space Flight Center (GSFC) over the past 6 years. It was originally known 
as the Low Cost Modular Spacecraft and Reference 1 provides an excellent 
description of its projected characteristics and capabilities. The brief 
discussion that follows is entirely based on Reference 1, and the illustra- 
tions of the MMS in Figures 2-2 through 2-5 are those of Reference 1. 

2.2.1 Multi -Mission Modular Spacecraft 

The MMS is an independent three-axis stabilized vehicle using three 
or four reaction wheels in conjunction with a stellar-inertial attitude 
reference system to provide accurate and stable spacecraft attitude con- 
trol. As the name implies, the MMS is composed of a set o f st andard mod- 
ules: ACS module. Communications and Data Handling (CDH) module and Power 
module. The propulsion module is generally optional and comes in two stan- 
dard versions known as SPS-I and SPS-II (L-D will require the SPS-I pro- 
pulsion module). The baseline MMS configuration contains the above men- 
tioned three subsystem modules supported by a module support structure, as 
shown in Figure 2-2. A transition adapter and a vehicle adapter complete 
the structural elements of the system. Mission unique and adapted items 
include the solar array, the solar array drive, the solar drive deployment 
mechanisms, a set of deployable TORS antennas and booms, and a mission 
adapter and payload (see Figures 2-3 and 2-4). The three spacecraft mod- 
ules are physically the same size. As shown in Figure 2-3, they measure 
18 inches deep, 48 inches high, and 48 inches wide. The ACS module is of 
particular interest in this study and a preliminary equipment layout of 
the module is shown in Figure 2-5. 
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PAYLOAD 





Figure 2-4. Multi -Miss ion Modular Spacecraft 
(Taken from Reference 1 ) 
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TORQUER BARS 





REACTION WHEELS 

REACTION WHEEL 
(REDUNDANT-OPTIONAL) 

INTERFACE ASSEMBLY 

REMOTE DECODER 
REMOTE MULTIPLEXER 

DIGITAL SUN SENSOR 


Figure 2-5. ACS Module 



Alignment between the attitude control subsystem and the payload ref- 
erence axis is maintained by the transition adapter during orbit opera- 
tions. The transition adapter has other unique features which allow a 
shuttle orbiter capture, docking, retrieval and servicing by a manipulator 
system. The transition adapter also supports the solar array launch re- 
straint and deployment mechanism and the solar array drive motors. The 
communication antennas are attached to the spacecraft to suit the mission. 
Optional attachment points include the area above the separation interface 
on the base adapter and also the transition adapter. The antennas and 
boom assembly can be folded into the recess formed by the adjacent edges 
of the modules. Total spacecraft weight without payload is 1235 pounds 
for the baseline configuration and 1590 pounds for the fully redundant con- 
figuration. This does not include the propulsion module. 

2.2.2 LANDSAT-D Configuration and Hass Properties 

Figure 2-6 shows a simplified block diagram of the L-D spacecraft sys- 
tem. On the left, the three standardized subsystem modules and their sup- 
porting structure is indicated. In the middle of Figure 2-6, the L-D pay- 
load is indicated. The payload instruments are supported by a mission 
unique forward structure and feed their output data to mission unique 
transmitters and antennas operating at K-band to handle the 120 megabit/ 
second TM data rates. Capability is planned for real time relay via the 
Tracking and Data Relay Satellite System (TDRSS) to the continental U.S. 
and for direct downlink to K-band equipped ground stations in view of the 
spacecraft. 

Figure 2-7 shows the overall L-D/MMS spacecraft in its on-orbit con- 
figuration with antennas and the solar array deployed. For this mission 
a one sided solar array is used to maintain clear fields of view for the 
passive radiation coolers used by both the TM and the MSS. For in flight 
replacement by the shuttle, the TM can be slid upward away from the nadir 
direction after disconnecting three attach fittings. The TDRSS compatible 
K-band antenna dish and transmitter are gimballed and deployed on a boom 
to allow the antenna more rotational freedom in tracking TDRS and to 
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Figure 2-6. Block Diagram of Assumed LANDSAT-D Spacecraft 
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minimize antenna radiation pattern interference with the solar array and 
other parts of the spacecraft. The direct access antenna for downlink 
K-band communication is not shown in Figure 2-7. 

The driven solar array is of the FRUSA type (Flexible Rolled-Up Solar 
Array developed by Hughes Aircraft Company) and it is canted 37.5 degrees 
relative to its axis of rotation (parallel to the spacecraft pitch axis) 
so that a coincal array envelope results. Because of the 9:30 a.m. orbit 
of L-D, the cant angle must be 37.5 degrees in order to maintain the plane 
of the solar array normal to the sun line when the array rotates about the 
spacecraft pitch axis at orbit rate (see also Figure 2-1). The array di- 
mensions are 6.33 ft x 27 ft (1.93 m x 8.23 m) of which the solar cell 
area comprises 5 ft x 27 ft (1.525 m x 8.23 m). The array weight is ap- 
proximately 50 pounds (22.7 kg) and the array efficiency after 3 years on- 

p 

orbit is about 7.5 watts/square-foot (80.73 watts/m ), resulting in a 
little over 1 Kw of generated electric power. 




The total spacecraft on-orbit weight is estimated at 
m = 3500 lb (1587 kg) 

and preliminary assessments of the principal moments of inertia are (Ref- 
erence 2) ' ( 


1500 slug-ft^ 

(2034 kg m^) 

1970 slug-ft^ 

(2671 kg m^) 

1750 slug-ft^ 

(2373 kg m^) 


2.3 Payload Description 
2.3.1 Thematic Mapper 

The Thematic Mapper is planned as the primary payload instrument for 
LANDSAT-D, to be launched early in 1981. 

The name "Thematic Mapper" derives from the intended use of the in- 
strument to produce land use theme maps which may be tailored to the data 
requirements of particular researchers, responsible government agencies, 
or private companies who can make use of remotely sensed data to better 
manage our earth resources. The Thematic Mapper will provide significant 
improvements in instrument performance compared to the current generation 
Multi -Spectral Scanner (MSS). The Thematic Mapper will provide coverage 
in at least seven spectral bands versus five in the MSS, 30 meter spatial 
resolution versus 80 for the MSS, better radiometer accuracy, and better 
geometric accuracy. Figure-2-8 shows a conceptual drawing of LANDSAT-D 
in operation. The spacecraft is in a circular, nearly polar orbit and 
provides complete global converage every 16 days. 

TRW Defense and Space Systems, teamed with Perkin-Elmer, and Hughes 
Aircraft Company (HAC) have submitted competing bids to NASA GSFC to de- 
velop and build a flight version of the TM. This competition has not yet 
been decided and for this reason the description of the TM provided here 
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Figure 2-8. Crosstrack Optical Scanning Provides Continuous 
Coverage of 185 km {100 nmi) Swath 


must be more of a general nature, omitting details. The Thematic Mapper 
optics include a scan mechanism which causes the line of sight of the in- 
strument to scan across track; with the TRW/P-E design each scan is linear 
and from West to East. The scanned image is relayed to a focal plane as- 
sembly where the light is split into at least six spectral bands (options 
may add a seventh and eighth band). Each 0.004 inch square detector in a 
line array has a 30 meter square footprint on the ground when projected 
back through the optics. The six spectral bands are specified as follows 
where band edge is the wavelength that results in 50% of the maximum re- 
sponse. 
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Band 

Number 

Lower Band 
Edge i-m) 

Upper Band 
Edge (ym) 

Maximum Edge 
Slope (ym) 

1 

.0.45 + 0.01 

0.52 + 0.01 

0.02 

2 

0.52 + 0.01 

0.60 + 0.01 

0.02 

3 

0.63 + 0.02 

0.69 + 0.01 

0.02 

4 

0.76 + 0.02 

0.90 + 0.01 

0.03 

5 

1.55 + 0.02 

1.75 + 0.02 

0.05 

6 

10.4 + 0.10 

12.5 + 0.10 

0.30 


In each spectral band, the moving image scene is swept past a line 
array of detectors (with the line array parallel to the direction of 
flight) by the scanning action crosstrack as outlined above. As the 
ground footprint of each detector moves 30 m across track, the detector 
output is sampled, converted to an 8 bit digital word, and ultimately 
multiplexed into a single 120 megabit/second data stream which is relayed 
to the wideband communications module. 

As indicated in Figure 2-8, each scan sweeps out a swath 960 meters 
(along track) by 185 kilometers (cross track) and the scan rate is chosen 
so the swaths just touch, providing continuous coverage along the space- 
craft ground track. 

The major difference between the TRW and MAC TM designs is the scan 
mechanism. The HAC design uses a reciprocating scan mirror while the TRW 
design uses a uni-directional scan wheel, referred to as the roof mirror 
scanning wheel, since it has roof mirrors spaced around its circumference. 
The roof mirror scanning wheel converts the fixed image from the telescope 
into an image which moves across the detectors. The concept i's illustrated 
in Figure 2-9. 

The scan mechanism of the TM is of some importance for the ACS studies 
because of potential attitude disturbances resulting from the moving parts. 
During normal on-orbit operations the TRW scan wheel rotates at a constant 
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Figure 2-9. Schematic of the Scanning Telescope Optical Path 


2 

Speed of 23.45 RPM and has a moment of inertia of 224 ounce-inch-seconds 
(1.58 kg-m ). No attitude disturbance occurs in steady state operation. 

The angular position of the HAC reciprocating mirror follows a triangular 
wave shape with rounded corners at a frequency of 7.14 Hz and 4^74 mradian 
(4.24 degrees) amplitude. The moment of inertia of the mirror about its 
pivot is 0.02614 kg-m^. 

The TM is nominally in a data taking mode only during those portions 
of the orbit when the L-D is flying over land masses which are in daylight. 
Data obtained over water (oceans) or a dark earth are of little practical 
value and consequently this type of data will be acquired only occasionally 
for experiment's sake. The fact that nominally no TM data need be trans- 
mitted to TORS over water or over a dark earth is an important ground rule 


2-15 



in this study. It is key to the proposed solar array control strategy end 
the proposed TDRS antenna gimbal configuration and mast length. The pro- 
posed configuration does allow, however, for data being taken at night and 
over water, but not under any arbitrary circumstances. This will be dis- 
cussed in more detail in Section 3.0, 

To give a feel for the size of the TM instrument the following ma'ss 
property and power data of the TRW design are provided. 

Thematic Mapper Physical Parameters (TRW Design) 

• Outline Dimensions: 1,12 x 1.94.x 0.96 meters 

(3.7 X 6.4 X 3.1 ft) 

• Overall* Weight: 223 kg (490 lb) 

• Power Consumption: 228 watts (operating mode) 

152 watts (night only) 


2.3.2 Multi -Spectral Scanner 

The secondary payload for L-D is the Multi -Spectral Scanner (MSS) as 
it will be flown on LANDSAT-C. The MSS is currently flying on LANDSAT's 
A and B, but the LANDSAT-C version will contain an additional (fifth) 
channel . 

The MSS is an object plane scanner that was developed by Hughes Air- 
craft Company to provide radiometric measurements and imagery of the 
ground scene in the following spectral bands: 

(1 ) 0.5 to 0.6 microns 

(2) 0.6 to 0.7 microns 

(3) 0.7 to 0.8 microns 

(4) 0.8 to 1.1 microns 

(5) 10.4 to 12.6 microns (LANDSAT's Cand^D only) 
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The uni -directional cross track scanning is provided by a high duty 
cycle (65%) rocking mirror at a frequency of 18.5 Hz over + 74 mradians 
(4.24 degrees). The total weight of the MSS is 125 pounds (56.7 kg) and 
the power consumption is about 60 watts. The MSS data rate is 15 megabits/ 
second. 

2.4 LANDSAT-D Pointing Performance Requirements 

Preliminary spacecraft pointing requirements are specified as follows: 

• Pointing Accuracy 

Roll and Pitch: 0.01 deg (lo) 

Yaw: 0.05 deg (la) 

f Pointing Stability 

Average Rate/Axis: < 10~ deg/sec (la) 

Stability over 20 Minutes: 

Roll and Pitch: +10.5 yrad (la) 

Yaw: ± 80 yT'ad (la) 

These requirements are referenced to the spacecraft coordinate reference 
axes which coincide with the optical cube on the star tracker. This means 
that the spacecraft to TM misalignments are not included in these perfor- 
mance specifications. 

The ultimate temporal registration accuracy desired (essentially re- 
peatability of the ground track after 16 days) for the TM is 3 meters 
(0.1 of a pixel) which corresponds to an attitude stability of about 4.2 
yrad over 16 days if no ephemeris error were incurred. Work on the 
autonomous on-board Thematic Mapper boresight registration is still being 
pursued and the study results will be presented in the final report. It 
is clear, however, that registration within 3 meters is a very stringent 
performance goal for LANDSAT-D. 
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3.0 PAYLOAD DATA RELAY SYSTEM 


This section addresses various system aspects of the TM data relay 
system. In particular, system configuration, control, and operational as- 
pects of the mast-deployed, gimballed Ku-band antenna, relaying TM data to 
TDRS, will be discussed. The antenna will often simply be referred to as 
the TDRS antenna (of L-D) to distinguish it from the undeployed and ungim- 
balled- direct access antenna which must also operate at Ku-band. This 
section will also discuss the interference of the solar array with the 
TDRS antenna beam and propose a solution to the problem. A digital comput- 
er program that computes on-orbit antenna motion profiles, determines per- 
cent solar array interference incurred, and determines whether the space- 
craft is over land or sea, over a lighted earth or a dark earth, is also 
presented. 

In order to decide whether the LANDSAT-D TDRS antenna can be pointed 
open loop by programmed command, or whether closed loop RF tracking, i.e., 
autotracking, is required, an antenna pointing error analysis was conducted 
first. This is the subject of the following section. 

3.1 TDRS Antenna Pointing Error Analysis 

Due to the constraint of limited TWT power, a large size L-D antenna 
dish (6 feet) is required in order to achieve sufficient EIRP for the re- 
turn link to the TDRS. This results in a narrow beamwidth of about 0.7 de- 
grees and necessitates very accurate antenna pointing in the order of 0.1 
degree (3a). 

This section discusses the antenna pointing accuracies available from 
two basic antenna tracking techniques: program pointing by ground or on- 
board stored commands (open loop) and autotrack (closed loop). Detail 
sources of pointing error are identified for both systems and the associ- 
ated errors are computed or estimated from accumulated data on current pro- 
grams and tracking systems hardware. 
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The results show that the proposed open loop pointing error (Rayleigh 
circular error probability) is 0.61 (3cr), which corresponds to an antenna 
gain loss (from peak gain) of 8.2 dB for receive, and 9.9 dB for transmit, 
and that the proposed closed loop pointing error (Rayleigh CEP) is 0.052 
degree (3o), which corresponds to approximately 0.1 dB antenna gain loss 
for both receive and transmit. 

The major error source for the open loop system comes from thermal 
distortion of the assumed 20 foot "Astro Mast" type of antenna boom. Since 
the closed loop system can track out the pointing error caused by antenna 
boom distortion, as well as many of the other-open loop pointing errors 
identified in the open loop system budget, a considerable reduction in 
pointing error is realized with the recommended autotrack antenna system. 

3.1.1 Program Pointing Performance 

Six general sources of error were found to contribute to the inaccu- 
racies involved in open loop pointing the L-D Ku-band antenna at the TDRS. 
These error sources, their breakdowns, and the error magnitudes are summa- 
rized in Table 3-1 . 

Spacecraft Position 

The errors in the target direction arise from uncertainties in the 
position of the TDRS and the L-D. The maximum error in the TDRS position 
is 61 m, assuming there is no correction immediately following each TDRS 
momentum dump (Reference 3). The TDRS orbit radius is 4.224 x 10^ km and 

3 

the L-D orbit radius is 7.083 x 10 km, so the minimum distance between 

4 

them is about 3.5 x 10 km. At this distance, a 61 m error is equivalent 
-4 

to a 10 degree 1 ine-of-sight error. An L-D position uncertainty of 
100 m (Reference 1) contributes 1.64 x 10”^ degrees error. These two ef- 
fects are negligible compared to the other errors shown in -Table 3-1. 

Spacecraft Attitude 

The spacecraft alignment error is dominated by the 0.01 .degree point- 
ing accuracy of the spacecraft and the relative error at the canister mount- 
ing surface due to thermal and other distortions of the spacecraft itself. 
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Table 3~1. Antenna Open Loop Pointing Error Sumnary 


Error Source 

3o Error (degree) 

1 . 

Spacecraft position 


1.92 X 10-^ 


TDRSS ephetneris, 61 m 

O 

X 

o 

1 



LF/0 ephemeri s , 100 m 

1.64 X 10"^ 


2. 

.Spacecraft attitude 


0.014 


Attitude control of ACS module (gyros, 
star tracker and RW control errors) 

0.010 



Spacecraft structural/thermal align- 
ment (ACS to canister mount) 

0.010 


3. 

Pointing computation and command 


0.023 


Computer roundoff and quantization 
(16 bits) 

0.005 



Gimbal angle encoder (14 bits) 

0.022 


4. 

Servo and gimbal led platform errors 


0.023 


Servo step size (0.03 degree) 

0.015 



Gimbal axes alignment (orthogonality 
and bearing runout) 

0.017 


5. 

Antenna mast errors {20' Astro-mast) 


0.520 


Canister mounting and alignment 
(thermal ) 

0.017 

\ 


Mast deployment/repeatability 




Bending 0.1“ (3o) 

Torsional 0.12“ {3o) 

0.129 



Thermal distortion 




Bending 0.05“ (3o) ^^p 

Torsional 0.5“ (3o) 

0.503 


6. 

Antenna errors (boresight relative to 
mounting platform) 


0.093 


Antenna attachment misalignment 

0.057 



Mechanical (thermal) boresight alignment 

0.028 



Electrical boresight alignment 

0.057 



Precomparator amplitude imbalance 

0.036 



Precomparator-postcomparator phase 
imbalance 

0.016 



Total rss error (3a) 


0.529 
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Pointing Computation and Command 

The computation and command error is dominated by the gimbal angle 
encoder error. The on-board data handling system uses 16 bits in double 
precision, but the payload drive encoder may well have fewer bits. The 
computational error, is found by letting the N bits cover a complete 
360 degree range 


_ 360 deq 

The computer roundoff and quantization error is found by letting N = 16 or 
-3 

e^g = 5.49 X 10 degrees. The gimbal encoder error is e^^ = 0.022 degree. 
The resultant error is 0.023 degree. 

The gimbal, servo, and platform error capabilities were taken from 
Reference 4 where it was determined that the combined gimbal and steady- 
state servo error is 1 arc-minute (0.017 degree) and the error due to gim- 
bal step size of 0.03 degree is 0,015 degree. The total error for this 
category is 0.023 degree. 

Antenna Mast Errors 

The main source of pointing errors is the Astro-mast. Basically, 
these errors relate the total error from the canister mounting surface to 
the gimbal mounts. The baseline mast is assumed to be 20 feet long* and 
9 inches in diameter when deployed. The mounting and alignment error in 
the table (0.017 degree) is the mast mounting error. The other errors are 
based on information obtained from Mr. Peter Preiswerk of Astro Research. 


*It will be shown in subsequent sections that this was a conservative 
assumption, i.e., the actual mast length is expected to be about 10 feet 
or less. The conclusions arrived at, namely that autotracking is required, 
is still believed to be valid, however. 
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■ The 3a numbers provided by Hr. Preiswerk assumed a standard 9-inch 
diameter, 20-foot long Astro-mast that is canister deployed. These num- 
bers are given in Table 3-1 under item 5. The bending and torsional er- 
rors are combined and represented by a Rayleigh circular error. The maxi- 
mum antenna pointing error due to the Astro-mast is 0.520 degree. 

Antenna Errors 

The antenna errors category includes all alignment errors from the 
antenna platform to the beam center. The 0.093 degree error assumes a 
graphite reflector and the error is compatible with that of a similar 5- 
foot dish currently being used on another TRW project. 

The resultant RSS errors for open loop pointing is 0.529 degrees/axis 
or a Rayleigh circular error of 0.61 degree (3a). Since an open loop 
pointing accuracy of about 0.1 degree (3a) is required for satisfactory 
communication with TORS, -it must be concluded that an autotrack system is 
required. 

3.1.2 Autotrack Performance 

A description of the recommended antenna autotrack system and acquisi- 
tion procedure is presented in a companion report on the L-D mission unique 
communication system, Reference 5. Autotrack performance is given there 
in terms of the RSS of tracking errors due to thermal noise, precomparator 
amplitude imbalance, precomparator and post-comparator phase imbalances, 
and drive step size. The autotrack errors are summarized in Table 3-2. 

The single-channel amplitude comparison autotrack system shown func- 
tionally in Figure 3-1 has been recommended over other tracking techniques 
such as interferometer, self (beam)-steering, mechanical scan, conical 
scan, etc., because of its low risk, low cost, minimum hardware complexity, 
and its good performance in tracking the Ku-band, spread spectrum RF signal 

In summary, the recommended autotrack design provides excellent track- 
ing performance. Practical limits of precomparator amplitude imbalance 
and precomparator and post-comparator phase imbalances do not cause 
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Table 3-2. Auto track Errors (3a) 


Error Source 

Error (deg) 

Notes 

Precomparator amplitude imbalance 

0.036 


Precomparator/ post-comparator 
phase imbalances 

0.018 


Drive step size 

0.015 

0.03° step size 

Servo loop dynamics 

0.015 

Estimated 

RSS 

0.0455 


Thermal noise 

0.0004 

(S/N)jp = 11 dB 

RSS total 

0.0455 

Per axis 

Circular error 

0.052 


Antenna pointing loss 



Receive 

0.06 dB 

3 dB beamwidth = 0.80° 

Transmit 

0.07 dB 

3 dB beamwidth = 0.73° 
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appreciable contributions to the autotrack error. The errors associated 
with motor step size and servo loop dynamics are also minimal. The track- 
ing error due to thermal noise is negligible at acquisition and autotrack 
receive signal levels. 

3.2 Antenna Gimbal Configuration 

Figure 3-2 shows a sketch of L-D in its normal on-orbit orientation. 
The figure mainly serves to define the body-fixed reference coordinate set, 

A. ^ A 

{Xfaj Zjjls and the earth half cone angle which equals 64.2 degrees 
at a spacecraft altitude of 705 km. It follows from the geometry indicated 
that the gimbal led Ku-band antenna must be able to track TORS in a cone of 
approximately + 120 degrees (exactly: 90° + 25.8° = 115.8 degrees) cen.- 

A 

tered about the -z^ axis, and that the maximum antenna dip angle below the 
spacecraft's local horizontal plane is 25.8 degrees. 

The Ku-band antenna target, namely TORS, follows a cone about the 
spacecraft’s y^^ axis which is constantly changing in size. This means that 
a two-gimbal arrangement should have the outer gimbal axis in or near the 
X|^-Z|^ plane so that frequent gimbal resetting is not necessary. The prob- 
lem with a two-axis (2~degree-of-freedom) gimbal lock problem exists and 
its effects must be carefully considered in the choice of the gimbal con- 
figuration and the degress of freedom. Four different gimbal configura- 
tions were considered. They are: 

(1) X over y* two-axis gimbal. 

(2) y over x two-axis gimbal. 

(3) Azimuth-elevation two-axis gimbal. 

(4) Three-axis gimbal using azimuth-elevation plus a third axis, 
cross elevation gimbal. 


*The notation means x rotation after y, i.e. , the outer gimbal axis is 
aligned with the S/C yj^ axis and the inner gimbal axis (at outer gimbal 
null) with the S/C Xj^ axi^s. 
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Figure 3-2. LANDSAT-D Normal On-Orbit Attitude and Geometry (not to scale) 


The X and y two-axis gimbal would need +180 degrees freedom about 
the y-axis and 0 degrees to 120 degrees in x. Since the line-of-sight to 

A 

TORS tends to come around the spacecraft's yj^ axis, this gimbaT configura- 
tion would operate often close to gimbal lock and high rates about the y 
gimbal axis would be required. For the same reason, it would also require 
frequent resetting of the outer gimbal because of the effective + 180 de- 
gree stops necessitated by the coax cables which are run across the gim- 
bals and do not permit a wind-up motion about y^. Since the gimbal lock 
tends to interfere with the L-D operational scenario, the x over y gimbal 
arrangement is discarded as a configuration candidate. 

The y over x gimbal would need 0-120 degree freedom in both gimbal 
axes. Gimbal lock can occur twice per orbit. Gimbal lock would occur when 
the TORS is either straight ahead along the spacecraft velocity vector 
(+Xj^) or when the TORS is exactly behind the spacecraft (-x^^ direction). 
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These situations would occur when TDRS is crossing the L-D orbit plane at 
the time L-D descends from the North Pole toward the equator, and again a 
little later when L-D is half way on its way to the South Pole. The y 
over X gimbal is a feasible configuration but not particularly attractive 
since the gimbal lock conditions interfere with normal L-D operations. 

The azimuth-elevation gimbal (with azimuth being rotation about the 
Zj^ axis = mast axis) has one well defined gimbal lock condition per TDRS: 
Whenever TDRS is in the spacecraft zenith, i.e., whenever TDRS is in the 
L-D orbit plane directly above the spacecraft (-Zj^ direction). Fast an- 
tenna rotation near gimbal-lock occurs then mainly about yaw (Zj^ axis) and 
is expected to cause minimal interaction with the spacecraft which could 
adversely affect payload pointing. The mechanical design of an azimuth- 
elevation gimbal is relatively uncomplicated and, moreover, a third gimbal, 
should it ever be required, can be added quite easily by putting a cross 
elevation hinge into the antenna mast. The azimuth gimbal about should 
ideally have a freedom of 360 degrees without mechanical stops, but be- 
cause of a cable transfer across the gimbal hinge, a cable wrap providing 
200 degrees of freedom with attendant overlapping + 200 degrees mechan- 
ical stops, will be required. The elevation gimbal freedom must be at 
least from -26 degrees to +90 degrees. The azimuth-elevation gimbal con- 
figuration looks attractive and will be further considered. 

A three-axis gimbal avoids gimbal-lock conditions but the gimbal con- 
figuration is mechanically more complex, bulkier, heavier, and costlier. 

At certain times two gimbal axes may be parallel and control laws and gim- 
bal servo cross-talk problems become more complex. It has been decided to 
avoid a three-axis gimbal system if satisfactory operation can be achieved 
with a two-axis gimbal. 

3.2.1 Gimbal Lock Problem 

The potential gimbal-lock problem for the two-axis gimbal configura- 
tions in general, and for the azimuth-elevation configuration’ in particu- 
lar, can be solved in one of two ways. Referring to the gimbal-lock point 
as the. keyhole, one can 
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(1) Show that the target (TORS) is never in the keyhold. 


(2) Place the keyhole so that the data which is lost, by the 
antenna not tracking TDRS, is not required. 

Option 1 can be viewed as: show that whenever one TDRS is in the key- • 

hole, it is acceptable to use the other one as the target. This is best - 
considered in terms of the target location. For an azimuth-elevation ar- 
rangement, gimbal-lock occurs when the target is in the L-D orbit plane. 

At that time the other TDRS is 130 degrees out of the orbit plane where it 
is hidden by the earth; thus there is no increase in coverage available by 
switching targets — in fact, since the other target would have to be 
tracked for the whole L-D orbit, there would be a larger reduction in 
coverage by switching targets. This exercise can be repeated for other 
two-axis gimbal configurations, but the result is the same. Since the 
outer gimbal axis must be near the plane, gimbal-lock will occur 

when a TDRS is near the L-D orbit plane. At that time, the other TDRS is 
50 degrees out of the L-D orbit plane and cannot be tracked throughout the 
L-D orbit. 

Option 2 is a feasible solution. Using the azimuth-elevation arrange- 
ment, aimbal-lock will occur when the L-D is directly beneath the TDRS, 
and so loss of the TDRS link will occur at 41 degrees and 171 degrees west 
longitude, i.e., beneath the station locations of the two TDRS's. These 
points are over the ocean where coverage is not essential. The loss of 
the TDRS link at 41 degrees west occurs close to the coast of Brazil, and 
at 171 degrees west there may be some Pacific islands of interest. Thus 
it is of some importance to precisely establish the size of the loss of 
coverage zone on the ground and to investigate if this zone could be 
shifted if required. 

The loss of coverage zone can be shifted by tilting the outer gimbal 
axis, which for the azimuth-elevation gimbal would mean to tilt the mast. 
This is demonstrated in Figure 3-3, where tilting the gimbal axis by the 
angle a, moves the loss of coverage zone by an amount y* The relationship 
between these angles is found from; 
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Figure 3*3. Moving the Loss of Coverage Zone by Tilting the 
Outer Gimbal Axis 


ana 


thus 


a = $ + Y 


h = r^^ sin Y = ' ''l ^ 

I 


a 


Y + tan 


( rj_ sin Y j 
ri ” ’^L '*') 



(3.2-1) 


(3.2-2) 


(3.2*3) 


where the approximation is for small angles. Using the values r, = 7074 km 

4 ^ 

and rj = 4.224 x 10 km, a « 1.2 y- For example, if the gimbal axis is 

tilted 6 degrees, the coverage hole moves about 5 degrees. 


The size of the loss of coverage zone can be estimated by assuming 
that normal operation of the antenna gimbal can be maintained until one 
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gets within + 2 degrees of gimbal-lock. Using Equation (3.2-3) with 
a = + 2 degrees, this becomes y ~ ±_ 1.664 degrees of the earth's sphere on 
the ground which equals + 185 km. Since a swath of the TM happens to be- 
exactly 185 km wide in the east-west direction, the loss of coverage zone 
in the worst situation is at most three swaths wide (555 km w 5 degrees on 
earth sphere), and for the best situation at least two swaths wide 
(370 km » 3.3 degrees). In the north-south direction, the loss of cover- 
age zone extends always + 185 km (+ 1.664 degrees) if gimbal-lock can be 
approached to within 2 degrees. This is illustrated in Figure 3-4. Re- 
ferring to appropriate maps showing the north-east coast of Brazil and the 
Pacific Ocean, respectively, it can be seen that the loss of coverage zone 
contains only water (the coast of Brazil is barely missed for the worst 
situation), with the exception of perhaps some of the Gilbert Islands in 
the Pacific. Thus, the antenna mast need not necessarily be tilted. If a 
tilt is considered, because normal tracking operation can only be main- 
tained to within, say, 5 degrees of gimbal-lock, then the mast could be 
tilted by 6 degrees to 12 degrees toward +Xj^ so that the now larger loss 
of coverage zone would move 6 degrees to 12 degrees north at the descend- 
ing node (daylight), where no nearby land masses exist. But this is ac- 
ceptable only if the loss of coverage due to gimbal-lock is acceptable at 
night: if the zone is shifted one way at the descending node, it is 

shifted in the opposite direction at the ascending node. It would be pos- 
sible to place the zone over the ocean all the time by tilting the gimbal 
axis 36 degrees from the axis in the plane. This moves the zone 
30 degrees north or south depending on the L-D direction at the time. 

3.2.2 Baseline Gimbal Configuration 

The azimuth-elevation gimbal configuration has been selected as the 
baseline gimbal for the L-D Ku-band (TORS) antenna. The occurrence of 
gimbal-lock is well defined for this gimbal configuration and it was shown 
that the resultant loss of coverage zone is always situated over water. 

The gimbal configuration is shown in its null position in Figure 3-5. 
Gimbal-lock occurs when the elevation angle g« = 90 degrees and the anten- 

A. Ax ^ Ca 

na views zenith. The coordinates {x, y, z} are identical to the body 
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Figure 3-4.“ Zone of Loss of TDRS'Line Due to Gimbal-Lock (Normal Operation 
' Assumed to Within +,2 degrees of Gimbal-Lock) 


3-14 


Zeni th 



Nadir 

Figure 3-5. Ku-Band Antenna in the Null Position 


•fixed set defined in Figure 3-2, i.e., in its null position the antenna 
faces in the direction of travel of the satellite. This is important in 
view of the cable wrap limits specified as + 200 degrees from the origin, 
since near gimbal-lock, the antenna must be able to rotate nearly 180 de- 
grees in azimuth from +x to -x without crossing a cable wrap limit. A 
single mechanical stop in azimuth is unacceptable since the antenna must 
be able to fine-track a TORS at any arbitrary azimuth angle. Thus, an 
overlapping cable wrap and attendant overlapping mechanical stops (to 
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prevent cable damage) are required. The peak gimbal rate occurs in azi- 
muth near gimbal-lock. It can be shown that if gimbal-lock is approached 
to within e degrees, then, for small e, the peak azimuth rate is given by 




9l 


(3,2- 


max 


where is the orbit rate of L-D. For e = 2 degrees, this yields 


g, = 1.74 deg/sec 
max 


The peak azimuth rate could be decreased if the elevation angle were al- 
lowed to "follow through" beyond 90 degrees, i.e., the elevation gimbal 
freedom were -30° 1 210 degrees. The software of the gimbal control- 

ler would also become more complicated, however, since gimbal steering is 
not unique anymore. 

Table 3-3 summarizes the baseline gimbal design requirements. A ser- 
vo loop controlled stepper motor has been assumed for the gimbal drive. 
Step size and other requirements are compatible with the antenna pointing 
requirements and the autotracking baseline system established in Section 
3.1. The slew rate and acceleration requirements assure that: 

• The azimuth can be reoriented by 180 degrees sufficiently 
fast (1 minute) to resume tracking TORS when it emerges again 
from the gimbal-lock zone, assumed at 2 degrees about the 
keyhole. (That is, the time to reorient the azimuth gimbal 
is limited by the time it takes L-D to cross the gimbal-lock 
zone. ) 

• The antenna can be slewed from one TORS to the other in less 
than 2 minutes. 
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• The system has sufficient acceleration capability for a 

spiral search pattern to find and lockup on TORS, in case the 
open loop pointing command did not bring the target within 
the acquisition field of view of the autotrack system. 


Table 3-3. Baseline Gimbal Configuration Design Requirements 


Angular Freedom; 


Azimuth 

+ 200° (overlapping stops) 

Elevation 

-26° to +90° (minimum) 
-30° to +210° (desirable) 

Maximum Tracking Rate: 

2.0 deg/sec 

Maximum Slew Rate: 

4.0 deg/sec 

Maximum Acceleration: 

2 

1.5 deg/sec 

Output Step Size: 

< 0.03 deg 

Readout Accuracy: 

Better than 0.2 deg 

Open Loop Pointing: 

Better than 0.25 deg 

3.3 Antenna Control System Functional Design 


This section presents a preliminary design for a control system to 
point the L-D Ku-Band antenna at a TDRS. There are two major control 
modes. Mode 1 is an open loop pointing mode in which the on-board com- 
puter generates gimbal angle commands based on the computed location of 
the target. In this mode the gimbal angle readouts provide the feedback. 
Mode 1 is used to: 
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(a) Slew from one TORS to the other. 

(b) Reset the azimuth gimbal angle when one of the + 200 degrees 
cable wrap stops are close to interf erring with nominal op- 
erati on . 

(c) Move the antenna through the zone of TORS link loss due to 
gimbal-lock. 

(d) Follow TORS through the zone of exclusion over the Indian 
Ocean where both TDRS's are hidden by the earth so that auto- 
track may resume as soon as a TORS comes again into the field 
of view. 

Operations (a) through (c) will nominally all be accomplished when the 
spacecraft is near the polar regions so that no data loss (i.e., no effec- 
tive communication loss) is incurred. The loss of coverage over the 
Indian Ocean, item (d) above, is not under the control of the L-D mission; 
it is a characteristic of the TDRSS as currently configured. 

If the open loop antenna pointing accuracy is worse than 0,657 de- 
grees, it may not be possible to establish the RF link for locking on to 
the TDRS. In this case the computer generates a spiral search pattern to 
bring the TDRS within the field of view of the autotrack system. The 
search pattern is generated in Mode la. 

Mode 2 is a closed loop tracking mode which does not use the on-board 
computer. The pointing error is sensed by the autotrack system, and the 
controller removes the measured error directly. 

3.3.1 Acquisition Procedure 

Acquisition procedures between the L-D an'd the TDRS are outlined in 
Table 3-4. The data was taken from the TDRSS "User Spacecraft Acquisition 
Procedures," P-805-1, June 1976, 6SFC. For L-D acquisition independent of 
the S-band TT&C link (except for initial antenna pointing programming and 
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Table 3-4. TDRSS User Spacecraft Acquisition Procedures 


3.1 KSA ACQUSITION SEQUENCE NO. 1 

3.1.1 An SSA return link is established with the user spacecraft using the normal 
SSA acquisition procedure. 

3.1.2 The user spacecraft radiates a minimum S-band EIRP of +13 dBw. 

3.1.3 TDRSS uses the S-band return link signal to aid antenna pointing. TDRSS 
radiates, within 5 seconds of user S-band turn-on and in the direction of 
the user spacecraft, a Ku-band +40 dBw (minimum) signal EIRP compatible 
with the forward link signal parameters defined in either Table 2-1 or 
Section 8. 2. 1.1 of S-805-1. The Ku-band carrier and PN clock shall include 
doppler compensation as scheduled. 

3.1.4 The user spacecraft acquires the TDRSS Ku-band signal, begins autotracking, 
and begins radiating a +30 dBw (minimum) signal EIRP at Ku-band in the 
direction of TDRSS. 

3.1.5 TDRSS then: 

a. Autotracks the user Ku-band signal within 5 seconds of user turn-on 
(P^cq ” ^-99 for user EIRP = +30 dBw). 

b. Begins radiating normal power mode signal EIRP or high power mode sig- 
nal EIRP as scheduled. 

c- Either establishes return link DG 1 service and/or DG 2 service, or 
Shuttle Ku-band service within 5 seconds. 

3.2 KSA ACQUISITION SEQUENCE NO. 2 

3.2.1 An SSA forward link is established with the user spacecraft using the nor- 
mal SSA acquisition procedure, 

3.2.2 The user spacecraft is commanded to point its Ku-band antenna at TDRSS and 
to transmit in the direction of TDRSS a minimum of +30 dBw EIRP at Ku-band 
compatible with the return link signal parameters defined in either Table 
2-2 or Section 8. 2. 1,2 of S-805-1, 

TDRSS searches for the Ku-band signal and begins autotracking within 10 
seconds of user turn-on (P,,^„ = 0.99 for user EIRP = +30 dBw). 

aCC) 

TDRSS radiates normal power or high power EIRP at Ku-band in the direction 
of the user spacecraft compatible with the forward link signal parameters 
defined in Table 2-1 or Section 9. 2. 1.1 of S-805-1, The Ku-band carrier 
and PN clock shall include doppler compensation as scheduled. 

3.2.5 The user spacecraft acquires the TDRSS Ku-band signal and begins 
autotracking. 

3.2.6 TDRSS either establishes return link DG 1 service and/or DG 2 service 
within 5 seconds. 


3.2.3 
' 3.2.4 
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other tracking considerations given below), KSA acquisition sequence 2 is 
recommended for the L-D mission. The acquisition procedure has also been 
discussed in Reference 5. 

L-D antenna tracking is accomplished as follows. Prior to coming . 
over the horizon, the L-D Ku-band steerable antenna is program-pointed at 
the TDRS within a factor of 1.8 (open loop pointing error) of the one- 
sided, half-power beamwidth (0.365 degree). This open loop uncertainty 
does not reduce the L-D EIRP below the minimum 30 dBW EIRP required for 
TDRS Ku-band antenna acquisition and autotrack. A margin of greater than 
18 dB exists. 

Upon acquiring and autotracking the L-D signal beacon, which consists 
of the normal composite 120 Mbps and 15 Mbps instrument signal, the TDRS 
begins transmitting a spread spectrum beacon signal to the L-D. The TDRS 
beacon signal will be within 10 dB of the L-D Ku-band high gain antenna 
peak. This signal level is above the threshold level which initializes 
autotrack. 

The autotrack design provides for the inclusion of a small acquisition 
horn (BW ~ 12 degrees) to provide a broadbeam Ku-band beacon for initial 
TDRS acquisition of the L-D. The acquisition horn is considered as a back- 
up acquisition aid in the event of temporary pointing anomolies (attitude 
control, boom distortion, etc.). Furthermore, a preprogram search routine 
in the communication and data handling module on-board computer provides a 
search pattern for L-D antenna acquisition and autotrack. This is also 
considered a backup mode. 

3.3.2 System Functional Description 

A functional block diagram of the controller is shown in Figure 3-6. 
The required electronics are shown functionally by the block diagram in 
Figure 3-6a; a separate compensator may be needed for Modes 1 and 2 due 
to the different requirements. The autotrack measures the antenna point- 
ing errors - and eg - in antenna-fixed axes and they need be resolved 
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Figure 3-6. Functional Block Diagram of Antenna Control System 
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NOTE: AND THE t^‘S MAY BE DIFFERENT 

FOR THE TWO MODES 


Figure 3-6a. Functional Block Diagram of Drive Electronics 






into gimbal coordinates. One of these axes is aligned with the inner gim- 
bal axis, so the measured pointing error is also the gimbal angle error, 
i.e., for the elevation gimbal. 


^2 = C 2 - 92 (3.3-1) 

where C 2 is the commanded gimbal angle and gg the actual (measured) gimbal 
angle. But, the pointing error about the outer axis (azimuth) differs 
from the gimbal angle error according to the target location: 

e-j = cos 92 (C.j “ g-j ) (3.3-2) 

This factor is removed as shown in Figure 3-6. 

In Mode 1, the computer must calculate the required gimbal angles to 
point the antenna at the target given the TORS location and the L-D loca- 
tion and attitude. The necessary equations will now be derived. 

Assuming that the L-D and TORS positions are known in the standard 
earth centered inertial (ECI) coordinate frame {1} (Zj axis through North 
Pole, Xj axis toward vernal equinox), denote the L-D and TORS position 
vectors by and respectively. Then the vector pointing from 

L-D to TORS is given in ECI coordinates by 

{p}j = (3.3-3) 

and in L-D body coordinates by 

{p)b = 

where is the direction cosine matrix of L-D specifying its attitude 
relative to ECI. Denote the normalized pointing vector to TDRS in L-D 
body coordinates by 
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(3.3-5) 


\ ='<pV11pI1 

Denote the direction cosine matrix that transforms from body fixed co- 
ordinates to antenna fixed coordinates. Then, performing two ordered ro~ 
tations g-j and 





cos 92 cos 
-sin 

sin 92 cos g-j 


cos 92 sin g^j 
cos 

sin 92 cos 


-sin 92 
0 

cos 92 


(3.3-6) 


Note that when g^ = 92 = 0, the antenna fixed axes coincide with the body 
fixed axes and the antenna boresight is aligned with the spacecraft x^^ 
axis. The gimbal angles g-| and 92 that align the antenna boresight with 
the normalized pointing vector Tj^, must then satisfy 


■ 1“ 


1 

• — 1 

cr 

X 

1 

0 

" \b \ ■ ^ab 

1 

<< 

_0_ 




(3,3-7) 


which leads to 


azimuth: 


S1 = 


and 


elevation: 92 « -sin”^ (T|^^) -30“ £ 92 90° 


(3.3-8) 


(3.3-9) 


When the antenna is pointed at the computed target location, Mode 2 
is initiated if the autotrack "sees" the TORS beacon, and Mode la 
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otherwise. In Mode la, the computer generates the search pattern and 
switches to Mode 2 when the autotrack "sees" the beacon. The search pat- 
tern assumed here is a spiral: 

sin 92 = K t sin o)t (3.3-10) 

^2 ~ ^2 ~ ^ cos tot (3.3-11) 

At the end of n full periods 


a-|{n) 


2TinK 

to 


sin 2 iTn 


0 


d2(n) 


2TrnK 

to 


cos 2 irn 


2TrnK 

to 


?tK 

Thus the distance between successive scan centers is d = — ^ . With the 

to 

0.8 degree FOV width, sufficient overlap could be achieved with d =0.6 
degree, or K/to = 0.0955 degree. This scan will cover a circle of radius 
r - 3 degrees in n = 5 cycles. 

The peak gimbal rate within that circle would occur at the end of the 
fifth cycle, where 


a-j(o) = [K sin tot + K tot cos tot] = 2 TtnK 

t = 

to 


(3.3-12) 


&o(o) = [K cos tot - K tot sin tot] = K 

to 


(3.3-13) 


are the rates at the end of the n 


th 


cycle. 


Further, from (3.3-10) 


9l - 



(3.3-14) 
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For acquisition with Qg 60 degrees » this becomes 

:i 0.5 g-j £ 1 deg/sec (3.3-15) 

• • 

where g-j and Q 2 are each limited to 2 degrees/second by the design require- 
ments. Combining (3.3-12) and (3.3-15) 

Kn £ 0.1592 deg 

Since this rate requirement is only necessary thru the fifth cycle, one 

chooses K = 0.03 degree/second. Since K/im ~ 0.0955 degree, this results 

in tu = 0.3T4 rad/second, and the period is 20 seconds/cycle. The peak 

gimbal acceleration can be shown to be 0.6 degree/second where the gimbal 

2 

design requirement is 1.5 degrees/second . In summary, in Mode la the 
computer must compute gimbal angles according to Equations (3.3-10) and 
(3.3-11) with K = 0.03 degree/second, and w = 0.314 rad/second. A computer 
cycle time of 256 msecond would allow the gimbal commands to be updated 
80 times per search cycle. 

3.4 Digital Program for Antenna Motion Profile and Solar Array Interference 

One of the requirements of the LANDSAT-D study was to develop a dig- 
ital computer program for studying solar array interference with the Ku- 
band antenna. Program LFO (Landsat Follow-On) satisfies this requirement 
and is additionally a more general program that can be used to study vari- 
ous other spacecraft/orbital phenomena. In particular, the main outputs 
of the program are: 

• Antenna gimbal angle motion profiles 

• Percent transmission to TORS (interference problem) 

• Spacecraft eclipsed or not eclipsed 

• Longitude and latitude of the subsatellite point 

• Subsatellite point over land or water 

t Subsatellite point over lighted earth or dark earth 
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The objective of this section is to discuss the underlying theory on 
which the program is based, provide a brief description of the program or- 
ganization and its use, and to present and discuss some of the results ob- 
tained. Program flow diagrams are given in Appendix A and a program list- 
ing (FORTRAN IV) is given in Appendix B. 

3.4.1 Program Description 

This program is primarily used for computations involving a spacecraft 
in an arbitrary circular orbit — specifically the LANDSAT-D (L-D) — and 
another spacecraft which is in an earth-synchronous orbit — specifically 
the TORS. The geometrical relationships are shown in Figure 3-7; this fig- 
ure also defines several of the variables and parameters used in the pro- 
gram. The earth centered, sun referenced, coordinate system, s, has the 

z axis pointing north along the polar axis and the sun lying in the y = 0 

plane generally in the +x^ direction; the x^ axis and y^ axis lie in the 
equatorial plane. The TDRS‘s are located by the angles y-j and measured 

relative to the x^ axis — nominally TDRSl is at 41 °W longitude and TDRS2 
is at 171°W longitude. The L-D orbit is defined by its inclination angle 
Yj and the location of its descending node yq- In "this orbit plane, the 
spacecraft is located by orbit angle yl tiieasured from the p"oint 90 degrees 
ahead of the equatorial plane. The orbit coordinate system is located at 
the spacecraft c.m. with the z^ axis pointing at the earth center and the 
Xq axis in the direction of spacecraft orbital motion. In the program the 
body coordinate system is assumed to be identical with the orbital frame 
0 — no pointing errors. The angle between the sun vector and the equa- 
torial plane is y i this varies over the range [-23.5 degrees, 23.5 degrees] 
depending on the time of year. 

The mainline program contains minimal computation. Instead, it pri- 
marily ties together the various subroutines; this allows the user to easily 
supply special purpose subroutines to meet the needs of a particular situ- 
ation. The flow charts of the total program are given in Appendix A. The 
subroutines DOT, CROSS, BASIS, CITOO, COTOA, UNIT and ROT are minor compu- 
tational subroutines used for various vector manipulations and coordinate 
system transformations; their specific functions are shown in Table 3-5. 
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Table 3-5. Minor Subroutines in LFO 


Subroutine 

Name 

Inputs 

Outputs 

Description 

DOT 

X, Y 

A 

A = I . Y 

CROSS 

A, B 

c 

C ^ A X B 

BASIS 

A 

B, C 

The unit vectors B and C form an 
orthonormal basis with A 

UNIT 

I 

Y 

Y = X/lIl 

ROT 

XT, A. N 

XO 

YO = [A]j^ XT , N is nth axis 

CITOO 


Y 

Transforms vector X in s coord, 
system to unit vector Y in 0 coord, 
system. 

. COTOA 

A, p 

B 

Transforms vector A in 0 coord, 
system to B in array^^tered 
coord, system (see Section 5). 


The subroutines INPUT and OUTPUT do just what their names imply; they are 
separate subroutines in order to keep the mainline program simple and 
straightforward. The remaining subroutines are the main computational ele- 
ments of the program; their functions are described briefly in Table 3-6. 
More complete descriptions of ARRAY and 6IMBAL are given in Section 3.4.3. 
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Table 3-6. The Main Computational Subroutines in LFO 


Subroutine 

Name 

Inputs 

Outputs 

Description 

^ UPDATE 


t, Y-|. Y2» Yl 

Increments variables between compu- 
tation times • 

I ARRAY 

T, S 

■R 

With antenna target at T and sun 
at S, R = Per cent of antenna pat- 
tern not blocked by array 

S ECLIPSE 

S, Rg 

I SUN, I DAY 

ISUN = y L-D in sun 

\ 0, L-D eclipsed 

( 1, subsatellite point 
IDAY = j lighted 

' 0, subsatellite point 
dark 

LAT 

l 

- 

Yr 

Yi> Yd* 
ALAT 1 

ALAT, ALONG 

Computes geographical latitude and 
longitude of subsatellite point. 

j GIMBAL 

T 

G,, Gj 

Given antenna target at T, computes 
the antenna gimbal angles 

1 LANDSEA 

i 

i 

ALAT, ALONG 

LAND 

LAND = 1 ?’ 

( 1) L-D over land 


The current version of the program assumes uniform step sizes between 
computation times, a simple rectangular array driven to follow the sun, 
circular orbits, and the- target (TORS) in a geosynchronous orbit. These 
assumptions are not necessary, and they could be changed by revising the 
various subroutines. For example, subroutine UPDATE could be changed to 
include the effects of an elliptical orbit or to make the intervals be- 
tween computation times vary — possibly dependent upon the program outputs. 
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Subroutine LAT contains an equation for converting true latitude to 
geographic latitude. The difference is shown in Figure 3-8 where (f)j is 
the true latitude calculated from the L~D position and c{>g is the geographic 
latitude which depends on the "local vertical". 



Figure 3-8. Illustration of the Difference Between True and 
Geographic Latitude 


In Reference 6, the two angles are related by the equation 


si n = 



(3.4-1) 
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where 


1.00336417 


b 

But this easily reduces to the more convenient form 

tan <})g = tan <f>y (3.4-2) 

which is used in the subroutine. This latitude correction changes the 
latitude by a maximum of 0.2 degree. The longitude computed in LAT is 
positive for western longitudes and negative for eastern longitudes. The 
subroutine LANDSEA uses the latitude and longitude to determine whether 
the point is over land or water. This subroutine was originally written 
for another program and has been incorporated into this program. The in- 
formation is stored in 42,261 bits, each representing a 111.2 x 111.2 km 
piece of the earth's surface — x P at the equator. 

3.4.2 Program LFO Abbreviated Users' Guide . 

The L-D user supplies the program input variables using a namelist 
file on TAPE3 called "IN". The variables supplied are listed in Table 3-7. 
The program supplies default values for most of the variables, as shown; 
the remaining variables must be supplied by the user. Some example runs 
are given later, but a sample input file is shown in Table 3-8. For the 
sample input shown, the output file — from TAPE4 — is shown in Table 3-9. 

In practice, the information on the output file would be plotted as in the 
examples given in Section 3.4.4. 

The input variables have the' units shown in Table 3-7. The various 
angles are supplied in degrees for convenience; the program converts the 
angles to radians for internal use, and then the output angles are recon- 
verted to degrees. 
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Table 3-7. Program LFO Input Variables 


Variable 

Name 

Symbol 

Description 

Default 

Value 

Alati 


Longitude of TORS #1 (Deg) 

+41.0 

ALPHAS 

a 

s 

Antenna FOV Skirt (Deg) 

2. 

■ ARRL 

£ 

Length of Solar Array (Ft) 

— 

ARRW 

u 

Width of Solar Array (Ft) 


ATMAST 

a 

Vector from Edge of Array to Base 
of MAST (Ft) 

0 

HH 

mm 

Solar Array Cant Angle (Deg) 

37.5 


■H 

Rate" ofl:hange of y-j (Deg/Hr) 

15. 


HI 

^Rate of Change of y^ (Deg/Hr) 

220.31 

DELTA! 


Time Increment (Hr) 


DIRflAST 


Vector Along MAST Toward Antenna 

(0,0,-l) 

FINAL! 

■H 

Final Time (Hr) 


6AMMA1 

H 

(Initial) Position of TDRS #1 (Deg) 

0. 

GAMMAD 


- Angle of Descending Node (Deg) 

37.5 

GAMI-lAI 

H 

L-D Orbit Inclination Angle (Deg) 

98.2 

GAMMAL 

(Initial) Position of L-D (Deg) 

0. 

GAMMAS 


Angle of Sun Below Equatorial Plane (Deg) 

0. 

HMAST 

h 

Antenna Mast Length (Ft) 


RDISH 


Radius of Ku-Band Antenna Dish (ft) 

3. 

REARTH 


Earth Half-Anqle from L-D (Deg) 

64,2 

RL 


L-D Orbit Radius (Km) 

7074. 

RT 

«T 

TDRS Orbit Radius (Km) 

4.22x10^ 


Table 3-8. A Sample LFO Input 


P5IM 

DT n AT *0 • 0 5 

FlNALT'l. 6 
H."‘A$7=t,34 

AF:RL*Z5 t6 
■ i£NO 
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Table 3-9. Output with Input Shown in Table 3-8 




T 

I 


J 



gp.EL 

^2 


9 |,Az 

92 

• ■ ■- 

TIME 

LP /0 

TD « S 1 

TRANS 

gimbal 

ANGLES 

TDRS 2 

TP ANS 

OiMBAL 

angles 

■ pay 

SUN 

LA 7 

lung . 

LAND 

6.000 

0.0 

0.0 

100.0 

- 37,2 

- 4,6 

"130 ,0 

ICO.O 

92.5 

- 17.4 

1 

1 

61,9 

- 11.5 

0 

.050 

11.0 

.7 

100.0 

- 39.1 

4,1 

- 129.2 

100.0 

90,2 

- 17,6 

1 

1 

76.4 

43,0 

1 

.100 

22.0 

1.5 

100,0 

- 42,2 

12.6 

- 126.5 

ICO.O 

87 iC 

-] 7.2 

1 

1 

66,7 

60,6 

1 

.150 

33,0 

2.2 

ICO.O 

- 46,8 

20,7 

- 127,7 

100.0 

85,7 

- 16.3 

1 

1 

56.2 

68,4 

1 

.200 

44.1 

3.0 

100.0 

- 53,1 

28.0 

- 127.0 

ICO.O 

03.9 

-14 ,e 

} 

1 

43.5 

73,1 

1 

.250 

55.1 

3.7 

100,0 

— 61 « 6 

34.1 

- 126.2 

100,0 

62.5 

- 13,0 

1 

1 

34 ,7 

76.6 

0 

.300 

66.1 

4,5 

100.0 

- 72,1 

38.3 

- 125.5 

100.0 

01.7 

- 10.9 

1 

1 

23.6 

79.4 

0 

,350 

77.1 

5.2 

100.0 

- 04,1 

40.2 

- 124,7 

• 100.0 

01.4 

- 8.7 

1 

1 

12,8 

81.9 

0 

.'tOO 

88.1 

6,0 

100.0 

- 96.2 

39.3 

- 124.0 

100,0 

81.7 

- 6.4 

1 

1 

1.9 

64.2 

0 

,< t 50 

99.1 

6.7 

100.0 

- 107.1 

35.9 

- 123.2 

100.0 

82.6 

- 4.1 

1 

1 

- 9.1 

06.6 

0 

.500 

110.2 

7.5 

100.0 

- 115.8 

30.7 

- 122,5 

100.0 

04.0 

- 2.1 

1 

1 

- 20,1 

89.0 

0 

.550 

121.2 

8,2 

100.0 

- 122.3 

24.2 

- 121.7 

100.0 

06,0 

-.5 

1 

1 

- 31.0 

91.7 

0 

.600 

132.2 

9 , 0 ' 

100.0 

- 126.9 

16.9 

- 121.0 

100,0 

00.3 

.7 

1 

1 

- 41,6 

94,9 

0 

.650 

143.2 

9,7 

100,0 

- 129.6 

9,3 

- 120,2 

100.0 

9 ] .0 

1.2 

1 

1 

- 52,6 

99,0 

0 

.700 

154.2 

10.5 

100.0 

- 131.4 

1.6 

- 119,5 

100,0 

93.7 

1.2 

1 

1 

- 63,2 

105.5 

0 

.750 

165.2 ■' 

11.2 ■ 

100.0 

- 131.8 

- 5.9 

- lie . 7 

100.0 

96.5 

. 4 ' 

1 

1 

- 73.3 

IIP . 2 

0 

.800 

176.2 

12.0 

100.0 

- 131.0 

- 13.2 

- 116.0 

100.0 

99.0 

-.9 

0 

1 

-01 ,C 

155 . e 

0 

.850 

167,3 

12.7 

91.4 

- 129,2 

- 20.0 

- 117.2 

100.0 

101.3 

- 2.9 

0 

1 

- 79,1 

- 137.0 

1 

.000 

196.3 

13.5 

- 20.0 

- 126.2 

- 25,6 

- 116.5 

100.0 

103.0 


0 ‘ 

1 

- 70.1 

- 1 X 1 . 4 ' ‘ 

1 

.950 

209,3 

14,2 

- 20.0 

- 122.0 

- 25,6 

- 115.7 

100,0 

104,2 


0 

0 

- 59.0 

- 101.5 

0 

1.000 

220.3 

15.0 

- 20.0 

- 116,5 

- 25,8 

- 115.0 

100.0 

104.0 


0 

0 

- 49.2 

- 96.0 

0 

1.050 

231 , 3 " 

15.7 

“ 20.0 

- 109.7 

- 25.8 

- 114.2 

100.0 

104.6 


0 

" 0 ■■ 

- 38. 4 

- 92.3 

'6 

1 1.100 

242.3 

16.5 

- 20.0 

- 102.0 

-2 5,8 

- 113.5 

ICO.O 

103,7 


0 

0 

- 27.5 

- 89.3 

0 

1.150 

253.4 

17.2 

- 20.0 

- 93,8 

- 25.0 

- 112.7 

100.0 

102.1 

- 20.8 

0 

0 

- 16*6 

- 86.7 

0 

1.200 

2 64.4 

18.0 

- 20.0 

- 05,7 

- 25.0 

- 112.0 

ICO.O 

99.7 

- 23,5 

c 

0 

- 5.6 

- 64 , 3 " 

0 

1.250 

275,4 

16.7 

- 20.0 

- 78.3 

- 25,8 

- 111*2 

100.0 

96.6 

- 25.6 

0 

0 

5.4 

- 02.0 

0 

1.300 

286,4 

19.5 

- 20.0 

- 72,1 

- 25.6 

- 110,5 

- 20.0 

92,9 

- 25.0 

0 

0 

16.3 

- 79,6 

1 

1.350 " 

297.4 ■ 

20,2 ■ 

- 20.0 

- 67.1 

- 25 . e 

- 109.7 

- 20.0 

8 8.6 ■ 

-2 5.8 

c “" 

■“0 

■ 27.3 

■■“-7 7. '6 

J 

1.^00 

308.4 

21.0 

- 20.0 

- 63,5 

- 25.8 

- 109,0 

- 20.0 

04,4 

- 25.0 

0 

0 

36.2 

- 74.0 

1 

1.450 

319.4 

21 .7 

48.3 

- 61.2 

- 22.7 

- 108.2 

- 20,0 

30.1 

- 25,0 

0 

0 

49.0 

- 70.3 

1 

1.500 

330,5 

22.5 

56.9 

- 59,9 

- 17,2 

- 107,5 

100.0 

76 . 0 ' 

- 25,3 

0 

1 

59,6 

- 64 . 9 ' 

.... 

1.550 

341.5 

23.2 

66.9 

- 59.6 

- 11*6 
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3.4,3 Two Key Subroutines 


The purpose of the subroutine ARRAY is to compute the fraction of the 
Ku-band antenna field-of-view pattern that is blocked by the solar array. 
The solar array is assumed to be a rectangular array of length sl and width 
w. It is also assumed that the array is mounted in such a way that it is 
canted to the spacecraft -y^^ axis and is driven to rotate about the y^ axis 
to track the sun. This is illustrated in Figure 3-9, where s is the sun 
vector, T is the target vector and 9 is the array cant angle. The dis- 
tance indicated as L in Figure 3-9 is the distance from the x-y plane to 
the top of the mast; this has been referred to as the "mast length" in pre- 
vious work. In the program, the variable HMAST is used instead to repre- 
sent the true mast length as measured from the spacecraft mounting point 
to the hinge. 

The field-of-view of the antenna is a cylindrical beam of radius R^ 
extending out from the antenna dish plus a 1 to 2 degree skirt all around 
the cylinder. The solar array interference is found by: 

(1) The solar array orientation is determined from the sun 
vector s. 

(2) A ray of the antenna beam at the center of the dish is found 
and used to determine the average distance from the dish to 
the plane of the array. 

(3) That distance is used to determine the FOV diameter at the 
plane of the array — this Includes the skirt. 

(4) A pair of unit vectors that form an orthonormal basis with 
T are used to construct a uniform distribution of points in 
the FOV. 

(5) A ray of the antenna beam through each of those points is ex- 
tended to the plane of the array. 

(6) The portion of those rays that intersect the array is the 
interference. 
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Figure 3-9. The Antenna and Solar Array in Orbit Coordinates 


Giver the sun unit vector, s = s , 

A y 

angle 6^ in Figure 3-9 is given by 



, the array rotation 


e 


s 



(3.4-3) 


where is the angle between the x axis and the projection of the sun 
vector onto the x~z plane. This means is zero when the spacecraft is 
over the North Pole — 8.2 degrees away, actually — at a solstice. If 
the array cant angle is then a vector in body coordinates 0 can be 
written in array coordinates A as 
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(3.4-4) 


'■a “ '"o = ■^A0'"o 

X y 

This assumes that the two coordinate systems have the same origin, which 
is an acceptable assumption in the case of the sun and TORS vectors. The 
coordinate system A is the one in which the solar array lies in the x^-y^ 
plane with its inner edgelying on the axis. 

In the program, the mast length h is called HMAST, the mast direction 
m is given by the vector DIRMAST and the antenna mast base is located by 
the vector ATMAST (?) relative to the origin of the coordinate system A — 
the above vectors are given in the 0 coordinate system. Thus, in array 
coordinates, the antenna hinge is at the point 

V ■= T^q {? + h m} (3.4>5) 


Since the antenna is pointed at the target, the center ray of the beam — 
which is aligned with the hinge — intersects the points 


r^ = v + a f (3.4-6) 

where a is a parameter. Thus the distance from the hinge to the plane of 
the array is given by 


ct 


0 



(3.4-7) 


because that is the value which makes (r^) = 0. This can be used to com- 
pute the FOV radius at the plane of the ar?ay: 




(3.4-8) 


where Rp is the radius of the dish, is the angle of the skirt and 
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d 


(3.4-9) 


LV 

tan '(^y 


where is the earth half-angle. 

Given — the target unit vector in array coordinates — and a suit- 
able unit vector a pair of orthogonal unit vectors in the plane P per- 
pendicular to are formed as 

®1 " ^ ^ ®o'^ (3.4-10) 


and 


e 


2 


■^A =1 


(3.4-11) 


These basis vectors are then used to form a set of 625 points uniformly 
distributed over the square in P that circumscribes the FOV circle of 
radius R . A particular point is given by 

c" = a^ e^ + a^ (3.4-12) 


where 


2i - n - 1 


R. 


i=l n 


(3.4-13) 


and 


"2 


2j - n - 1 
n 


R. 


j=l,...,n 


(3.4-14) 


where n = 25. 

A 

The ray parallel to T. through -one of the points is given by 
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r 


Fq + c 


(3.4-15) 


V + c + a 


If a is chosen so that (F)^ = 0, then r is the intersection of the ray with 
the plane of the solar array. This point of intersection is compared with 
the array boundaries to determine whether or not it intersects the array. 
The subroutine output is the percent transmission of the antenna which is 


R = 



(3.4-16) 


where M-| is the number of points lying in the circle and M 2 is the number 
of those points that emanate rays which intersect the solar array. 

Subroutine Gimbal 

The purpose of this subroutine is to compute the gimbal angles re- 
quired to point the Ku-band antenna at the target. In the subroutine, it 
is assumed that there are two gimbals with the inner gimbal angle being 
the elevation angle and the outer gimbal angle being the azimuth angle. 

In the normal or null position, the antenna is pointed along the roll axis 
and the elevation axis is the spacecraft pitch axis — see Figure 3-10. 
Thus for a normalized target vector T, the required gimbal angles g.j and 
92 are found from 


y ' z 

which is rearranged to obtain 



■fx' 


~cos g^ cos 

92 

T = 

Ty 

= 

sin g^ cos 

92 


Tz 


-sin §2 



1 

0 

0 


(3.4-17) 


(3.4-18) 
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Solving this for the gimbal angles: 


g^ = tan*"^ 

= sin"'' [-y (3.4-20) 

The gimbal equations for several two-gimbal configurations are given in 
Table 3-10. As discussed earlier, a third gimbal would be used only to 
overcome a severe gimbal-lock problem. Because, there are several ways to 
utilize the third gimbal, it is not convenient to include the three-gimbal 
configuration in Table 3-10. 




(3.4-19) 


Table 3-10. Antenna Gimbal Equations 


Outer 

Gimbal 

Axis 

Inner 

Gimbal 

Axis 

Nominal 

Antenna 

Direction 

Outer 

Gimbal 

Angle, 

9] 

Inner 

Gimbal 

Angle, 

92 

z 

A 

y 

[1 0 of 

tan"'*[Ty/T^] 

sin'’ [y] 

z 

y 

[0 0 -if 

tan-’EyT^] 

cos'’ [-T^] 

z 

X 

[0 0 -if 

tan"''[(-T^)/y 

cos'’ [y] 

A 

y 

X 

[0 0 -if 

tan"’[(-T^)/(-y] 

sin'^ [Ty] 

X 

y 

[0 0 -if 

tan"”* [Ty/(-T^)] 

sin"^ [-Tj^] 
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3.4.4 Exampi es 

Although the cases in this section are examples of the use of Program 
LFO, they are also meaningful cases to the L-D study. In all cases the 
antenna mast length above the attachment point has been chosen as 6.34 feet. 
In the first example, the L-D passes through the "zone of exclusion" — an 
area over the Indian Ocean where the earth blocks transmission to both 
TDRS’s simultaneously. The input file is shown in Table S-ll, and the re- 
sults of the run are shown in Figures 3-11 through 3-15. Figure 3-11 shows 
the percent transmission for the two potential targets. Note that toward 
the end of the first orbit shown, the earth blocks both targets — this is 
the "zone of exclusion" (ZOE). Figure 3-12 indicates when the L-D is over 
a lighted earth, when it is eclipsed,. and when it is over land. Note that 
the earth blockage of both TDRS's occurs over a dark earth in this case, 
but it could just as easily occur over a lighted earth. The solar array 
interference is much worse over this dark earth than over a lighted earth 
because the array follows the hidden sun. For the few orbits shown, TDRS2 
would be chosen as the target because of the 100 percent transmission over 
the lighted earth. 


Table 3-11. Input for First Example Showing 
the Zone of Exclusion 


PS IK 

DELTAT-0,02 
F1NALT*7 . 

HMAST=6.34 

ARRW«5 . 
ARRL=*25*6 

GAMKA1*jBj9 • 

SEND 
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.TRANSMISSION TO TORS 2 TRANSMISSION TO TORS 1 

-50.0 0.0 50.0 100.0 -50.0 0.0 50.0 




LF/0 ORBIT RNSLE ( DEG ) 


Figure 3-11. Zone of Exclusion. at NigRt' 


LP/0 OVER LRND LIGHTED LE/0 LIGHTED ERRTH 

-0.1 0,5 1.1 -0.1 0.5 1.1 -0.1 0.5 





Figure 3-12. Zone of Exclusion at Night 



ELEVRTION FOR TORS 1 ( DEG ) RZIMUTH FOR TDR3 1 t DEG ) 

-30.D 10. 0 SO.O SO.O -180.0 -60.0 60.0 160.0 




LF/0 ORBIT RNGLE ( DEG ) 


Figure 3-13. jone of ExcO usion at Night 







ELEVATION FOR TORS 2 ( DEO ) AZIMUTH FOR TORS 2 ( DEG ) 

r30.0 10.0 50.0 90.0 -180.0 -60.0 60.0 ' 180.0 




0.0 360.0 720.0 


1080.0 H40.0 1800.0 


LP/0 ORBIT RNSLE ( DEG ) 


Figure 3-14. Zone -of Exclusion at Night 



ELEVATION FOR TORS 2 ( DEG J ELEVATION FOR TORS 1 ( DEG ) 

-30.0 10.0 50.0 90.0 -30.0 10.0 50.0 90.0 



Azimuth for TDRS1 (Deg) 



Figure 3-15. Zone of txclusion at Night 
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Figures 3-13 and 3-14 show the required antenna gimbal angles. The 
elevation angle stops at -25.8 degrees because that is the edge of the 
earth: when the TORS position is at an elevation beyond -25.8 degrees, it 
is hidden behind the earth. Note in Figure 3-14 that TDRS2 passes through 
the L-D orbit plane, requiring that the azimuth gimbal be reset. Also, 
the elevation angle is very nearly 90 degrees — gimbal-lock occurs when 
the elevation angle is 90 degrees. Finally, Figure 3-15 shows the .eleva- 
tion angle plotted versus the azimuth angle. 

The results of the second example are shown in Figures 3-16 through 
3-20', the input file is in Table 3-12. This run shows the gimbal-lock 
case as indicated in Figure 3-20, According to Figure 3-16, the target 
should be TDRS2 during these four orbits. When the gimbal-lock occurs — 
for a L-D orbit angle of about 800 degrees as seen in Figure 3-19 — the 
earth blocks TDRSl , so TDRS2 is still the proper target. 

The third example shows the L-D passing through the ZOE over a lighted 
earth. The input file is in Table 3-13, and the results are shown in Fig- 
ures 3-21 through 3-25. Note that the solar array interference is much 
less than it was in the first example, this is because the solar array is 
rotated away from the antenna when the sun is overhead. 

A fourth example consists of a 24 hour run (15 orbits) at vernal 
equinox conditions. Of main interest here is a summary of the percent of 
lighted land encountered. The program was initialized so that the first 
orbit contained the longest continuous land swath extending from the 
northern rim of Russia, across the middle east, to the s.outhern tip of 
Africa. The lighted land data is important for the spacecraft power budget 
since nominally the TM needs be in the active data taking mode only when 
the spacecraft is over lighted land. The results are summarized in Figure 
3-26 which gives percent lighted land per half and full orbit and percent 
land per orbit. For this particular 24 hour run (repetition occurs only 
after 16 daysl) the L-D was over lighted land only 16,56% of the 24 hour 
period. The spacecraft was over land, lighted or dark, 34.23% of the 24 
hour period. 
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Table 3-12. Input for Second Example Showing 
Gimbal Lock 


PSIN 

DELT AT-O .02 
FINAL!*?. 

HMAST*6.34 

ARRW*5. 

ARRL*a5.6 

GAH^A1»36. 

SEND 


Table 3-13. Input for Third Example- Showing 
the ZOE Over a Lighted Earth 


PSIN 

DELI AT*0.02 
FINAL!*?. 

HMA$!*6.3A 

ARRW*5 • 
ARRL-25.6 
GANMA1»-159._. 
SEND 
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0.0 360.0 720-0 1080.0 H40.0 


LP/0 ORBIT ANGLE ( DEG ) 


1800.0 


Figure 3-16. Giinbal Lock 
3-50 


LF/0 OVER LRND LIGHTED LF/0 LIGHTED ERRTH 

“0.1 0.5 1.1 -0.1 0.5 I.l , -0-1 ' 0-5 



0.0 360.0 720.0 1080. 0 1440.0 1800.0 

LP/0 ORBIT RNQLE: ( DBG ) 


Figure 3-17. Simtial lock 





ELEVRTION FOR TORS 1 ( DEG ) RZIMUTH FOR TORS 1 ( OEG ) 

-30.0 ID.O 50.0 90.0 -180.0 -60.0 60.0 180.0 




LF/0 ORBIT RNGLE ( DEG ) 


Figure 3-18. Gimbal Lock 










' 0.0 360.0 720.0 1080.0 H40.0 1800.0 

LF/0 ORBIT RNGLE ( OEG ] 


Figure 3-19. Gimbal Lock 
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ELEVRTION FOR TORS 2 t DEG ) ELEVRTION FDR TORS I ( DEG ) 

-30.0 10.0 50.0 90.0 -30.0 10.0 50.0 90.0 



Azimuth for TDRSl (Deg) 



Figure 3-20. Gimbal Lock 
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TRANSMISSION TO TORS 2 TRANSMISSION T 



Lf/0 ORBIT RNGLE ( DEG ) 


Figure 3-21. Zone of Exclusion Over Daylight 
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LF/0 OVER LFIND LIGHTED LF/0 LIGHTED ERRTH 

0.5 1.1 -0.1 0.5 1.1 -O.l O.S 





LP/0 ORBIT ANGLE ( DEG ) 


Figure 3-22. of Exclusion Over Daylight 


ELEVATION FOR TORS 1 ( DEG ) AZIMUTH FOR TORS 1 t DEG )' 

-30.0 10.0 50.0 90.0 -180.0 -GO.O 60.0 180.0 




LP/O ORBIT RNGLE ( DEG ) 


Figure 3-23. Zone of Exclusion Over Daylight 
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ELEVRTION 

-30.0 



0,0 360.0 720.0 1080.0 H'lO.O 1800.0 

LF/0 ORBIT "angle ( DEG ) 


Figure 3-24. 7one of Exclusion Over Daylight. 
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Azimuth for TDRSl (Deg) 


o 



FiZIMUTH FOR TORS 2 ( Dl6 ) 


Figure 3-25. .Zone of dust on Over Daylight 
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Percent 

Lighted Land 
per Half Orbit 


Mean=16.56 

0 


Percent 

Lighted Land 
per Full Orbit 


Mean=16.56 




Percent 

Dark Land 
per Half Orbit' 




Percent 

Dark Land 
per Full Orbit 

Mean=34.23 


0 

0 1 








2/ ; ' r= 6 .|^7^ 8 D Ll 9 j : , "lb ^ U 

' ' - ' ' •’ ' !*, Number. of Orbits 


Figure 3-26. 24-Hour Run: Percent Liohted Land 


Note that in the percent per half orbit plots in Figure 3-26, the 
lighted or dark land does not always split exactly along the half orbit 
time line. The reason for this is the 98,2 degree orbit inclination and 
the half orbit definition as + 90 orbital degrees from crossing the equa- 
tor. On the light side of the orbit (descending path), for example, the 
spacecraft will reach dark land (the Antarctic) near its "southern" most 
point (90 orbital degrees after crossing the equator). This can be seen 
best on the "Percent Dark Land per Half Orbit" graph in Figure 3-26. 

3.5 Solar Array Control to Minimize Interference with Antennta Beam 

The purpose of this section is to support the recommendation of a 
solar array configuration and array control strategy such that optical in- 
terference with the Ku-band antenna beam is eliminated. This is to be ac- 
complished under the constraint of keeping the antenna mast as short as 
possible (mainly for spacecraft for structural reasons), since a suffi- 
ciently long mast permitting the antenna a field of view above all ob- 
stacles presents itself as an obvious solution to the problem. 

The L-D solar array can be stationary relative to the spacecraft, or 
it can be driven to rotate about the yj^ axis and follow the sun. The 
former design has the advantage of not requiring a drive mechanism, but 
it also has the potential disadvantage of necessitating a longer Ku-band , 
antenna mast. There are many options available in designing the solar 
array. Two assumptions were made to eliminate some of the unnecessary 
opti ons : 

(1) It is assumed that interference with TORS transmission over 
earth night is acceptable. 

(2) It Is assumed that the array shape will be rectangular. 

The spacecraft with the Ku-band antenna is shown in Figure 3-27. The 
antenna is a circular dish with a 0.91 meter (3 foot) radius mounted on 
and rotating about the end of the mast. The mast is mounted parallel to 
the axis, as shown, but generally lies near the edge of the spacecraft 
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Figure 3-27. Ku-Band Antenna Mounted on Spacecraft 



body rather than at the c.m. The mast length, L,‘ is measured from the 

A 

hinge point of the antenna; but note that the actual 
mast could be made shorter than this by mounting it on the top of the 
spacecraft instead of in the plane. The antenna field-of-view pat- 

tern is the 6 foot diameter cylinder extending out from the dish plus a 
1 to 2 degree skirt — the skirt is assumed to be 1.2 degrees in Figure 
3-27. The antenna is only allowed to dip 25.8 degrees below the space- 
craft horizontal, because that is where its line-of-sight grazes the earth — 
64.2 degree half angle. 


The region which does not interfere with the antenna FOV is the 63 de- 
gree half angle cone centered on the mast with its apex at the point where 
the edge of the dish touches the masf. This apex of the cone is then given 
by A = (0.91 m)/(sin 64.2°) = 1.01 meter (3.33 feet) below the hinge point. 
Define a new coordinate system by shifting the body coordinate system so 
that the origin coincides with the base of the mast. Now, if the 3-tuple 
(x, y, z) is used to specify a point in this coordinate system, the set of 
points lying within the cone satisfies the inequality 


+ y 


- 

\sin e^; 


+ z cot 0 


D 


(3.5-1) 


where is the earth half angle and is the maximum angle by which the 
antenna beam dips below the horizontal. 

3.5.1 The Driven Array 

The driven array will be a rectangular array canted 37.5 degrees to 

A A 

the -y^ axis and rotating about the yj^ axis so that its face remains nor- 
mal to the sun line. When the spacecraft is In the ecliptic plane over 
earth day, this array is in the position shown in Figure 3-28, which is 
the position of least interference with the Ku-band antenna. As the L-D 
rotates out of the ecliptic, the array rotates about the yj^ axis and moves 
toward the cone where it may interfere with the antenna. That interference 
is avoided by mounting the antenna at the end of a mast, but that mast 
should be as short as possible. 
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,ntenna 



Figure 3-28. The Spacecraft with the Driven Array 


The mast length can be minimized by proper choice of the shape of 
the rectangular array. The mast can also be kept short by keeping the 
array close to the position shown in Figure 3-28 while the antenna is 
transmitting. This can be accomplished by letting the array follow the 
sun whenever the spacecraft is over earth night but not eclipsed — inter- 
ference is acceptable -- and whenever the spacecraft is within some prede- 
termined angle, 9|^, of the ecliptic. When the spacecraft is in the 
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remaining region, the array will be stationary at the position which would 
have it aligned with the sun at 6j_ from the ecliptic (see Figure 3-29) » 

This approach reduces the average solar array power over an orbit. 

That average power is computed using 

-360 

^avg ^ OT J (3.5-2) 

0 


where 6 is the orbit position of the L-D and P{e)' Is the normal l 2 ed in- 
stantaneous power measured as the projected array area normal to the sun 
1-ihe divided by the nominal area 11.89 m^ (128 foot^). The L-D is eclipsed 
by the earth when within + 56.73 degrees of the ecliptic over the earth 
night, so if the array follows the sun throughout the orbit, 

{ 0, in eclipse region 

(3.5-3) 

1, outside eclipse region 

Thus, 

.213.27 V 

^0 = sfe I / N® = 0.6848 (3.5-4) 

*-33.27 ‘ 


for this ideal case, where the zero angle is taken to be the point of the 
orbit where the L-D is 90 degrees of revolution ahead of the daylight 
crossing of the ecliptic. This number, P^, represents the required average 
normalized power per orbit. 

Now, if the array is stopped when the L-D is more than e, from the 

2 ^2 

ecliptic, the average power per orbit with a 11.89 m (128 foot ) array 
would be given by: 
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0 + / 
ir J 

- 0 

cos (90 - e 

J 

e, + 33,27- 

sin (90 - 8, ) 


The efficiency is plotted versus e, in Figure 3-30. 

100 4 


Figure 3-30. Solar Array Efficiency versus Stop Angle from the Eclipti 


In, summary, the approach, wherein the array does not follow the sun 
throughout the orbit, reduces the average solar array power over an orbit. 

The overall reduction is on the order of 4% and relatively insignificant. 

It can also easily be compensated by a correspondingly small increase in 
array area. 

The minimum mast length can be determined by comparing the location 
of the worst-case point of the array with the cone represented by Equation 
(3.5-1). The worst case point on the array is some point on the "top" 
edge when the array is rotated the maximum amount toward the antenna — 6|^ 
degrees from the lowest point as shown in Figure 3-28. It is obvious that 
because the cone is convex the worst case point is at one end or the other 
of the top edge. Let the vector representing one of these points be "c. 

Then, from Equation (3.5-1), the minimum mast length so that the point F 
does not lie in the FOV of the antenna is given by 

^MIN " ■*'^y sin ~ (3.5-7) 

Let the array have length £, measured along the shaft, and width u. 

Now, the unit vector along the array shaft is 



■ 0" 


- sin sin e," 

£ = [e^ " 90]^ le^]^ 

-1 


c L 

’ cos * 


0 


sin cos e, 

c L_ 


(3.5-8) 


where is the array cant angle — nominally 37.5 degrees -- and the ro- 
tation [a]^ means a rotation of the coordinate system through the angle a 
about the v axis. The unit vector along the end of the array is 


to = 


- cos 

0 

- sin 6j^ 


(3.5-9) 
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Using these unit vectors, the vector representing an arbitrary point of 
the top edge of the array is 

c “ aA5.+ -^cj + a {3.5-10) 

a Z 

Where a^ is the vector representing the close end of the array shaft and 
a e [0,1] determines the particular point of the top edge — the two points 
of interest are and c’-j. The vector elements a^ and a^ are assumed to 
be 0, and Sy is nominally 0 but may be as much as -0.91 (3 feet). Equa- 
tions (3.5-6), (3.5-7) and (3.5-10) were used to create Table 3-14 showing 
the baseline design and several variations. 

It is evident by comparing the first and last entries of Table 3-14 
that the baseline design is much better than having the array follow the 
sun throughout the orbit. Note also that for a shallow array rotation 
angle such as = 45 degrees the FRUSA* array is nearly as good as the 
best rectangular array. Further, the table shows that only a slight in- 
crease in array area is needed to bring the efficiency up to 100%. Recall- 
ing that the actual mast length could be smaller than that shown by mount- 
ing it higher up on the spacecraft, it is concluded that a normal boom 
rather than an Astro-mast type boom is sufficient for mounting the antenna. 

This baseline design is recommended as the rotating array design. 

The array rotation scenario for an orbit is as follows: as the L-D comes 

out of eclipse, the array tracks the sun until the northern terminator is 
reached — this is within 23.5 degrees of the North Pole depending on the 
season; during this period the array could be in the antenna FOV, but there 
is no Thematic Mapper data to transmit. At the northern terminator, the 
array "jumps" ahead approximately** 45 degrees and stops until the L-D 
travels a compensating 45 degrees. At this time, the array is again driven 
to follow the sun through 90 degrees of L-D travel where it again stops 


*6.33 foot wide Flexible, Roll Up Solar Array with 5 foot wide cell region. 
**Incli nation i = 98.2° + 90°, 
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Table 3-14. Various Driven Array Configurations 


0 '^ 
a fet 


Mast-Array 

Separation 

(Fig. 2-2} 

Array 

Cant Angle, 
6(, (deg) 

Angle from 
Ecliptic within 
which the Array 
tracks the Sun 
0L ('J09) 

Maximum 
Antenna 
Dip Angle, 
00 (deg) 

Solar Array 
Actual 
Surface 
Area, 

Efficiency 

{%) 

Minimum Hast 
Length for 
FRUSA (m) 

(1.52 m = 5 ft 
wide cell reg.) 

Optimal 

Array 

Width 

(m) 

Minimum 

Hast 

Length 

(m) 

0 >. 

37.5 

45 

27 

12.32 

100 

2.11 

0.90 • 

1.82 

0 

37.5 

45 

.7 

11.89 

96.37 

2.11 

0.89 

1.80 

• 0 

37.5 

40 

27 

12.49 

100 

2.01 

7**“* 

— 

0 

37.5 

50 

27 

12.20 

100 

2.55 

i.4i 

i.35 

0 

30 

45 

27 

12.32 

99.14 

2.84 

1.85 

2.63 

0 

30 . 

45 

27 

■B 

100 

2.85 

1.85 

2.64 

0 

37.5 

45 

28 

12.32 

100 

2.28 

1.16 

2.06 

-0.61 

37.5 

45 

27 

12.32 

100 

2.38 

0.91 

2.09 

0 

37.5 

60 

27 

11.89 

98.91 

3.41 

2.04 

3.11 

0 

37.5 

70 

27 

11.89 

99.68 

4.29 

2.56 

3.70 

0 

37.5 

90 

27 

11.89 

100 

5.96 

2.75 

4.50 











































untn the L-D reaches the southern terminator at which point the array 
catches up. Finally, the array tracks the sun into eclipse and the pro- 
cess is repeated for the next orbit. This design insures that the array 
will not interfere with the antenna during transmission periods, but on 
many orbits the array can follow the sun throughout without interference 
due to the favorable TORS locations. This subject can be studied more 
thoroughly using the computer program presented in Section 3.4. Note also 
that while the L-D is eclipsed by the earth, the array position is imma- 
terial; thus it could be moved out of the way of the antenna to allow data 
transmission at night. 

3.5.2 The Fixed Array 

If the solar array is to be fixed relative to the spacecraft so that 
it can not follow the sun, it should always have the same area normal to 
the sun line. Since the sun cones about the vehicle pitch axis, the array 
should be mounted normal to this axis as shown in Figure 3-31. In this 
configuration, the solar array position can be moved along the axis. 
Thus, the only fixed distance to be considered is the distance from the 
lower edge of the array to the antenna pivot point; this is the dimension 
h in Figure 3-31. The antenna .position is assumed to be directly over 
the array in the baseline design. 


The largest rectangular array that can be enclosed by the triangle as 
shown has height a and base b where 


and 


a 


1 

2 



0.91' \ 
sin 0 ^/ 


b 


_J 

tan 0 JJ 



0.91 \ 
sin 8^/ 


(3,5-11) 


(3.5-12) 
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-ANTENNA 



Figure 3-31. Location of Fixed Array on Spacecraft 



Then the area is given by 


'' = “» = (h- (3.6-13) 

For this design, the sun strikes the array at an angle of 52.5 degrees, so 
the area must be 

A = {11,89 m^)/(cos 52.5°) = 19.53 (210.26 ft^) (3.5-14) 

and Equation (3.5-13) results in = 5.47 m (17.95 feet).. The lower 
edge of the array can not be much more than 0.61 m (2 feet) below the 
.Xb“yb plane without interfering with the direct access antenna. Thus the 
mast length, L, as defined earlier would have to have a length at least 
L = 4.88 m (16 feet) . 

This required mast length would be increased still further if the 
antenna were laterally offset from the array. Also, any change in the sun 
incidence angle would cause. a change in the array size and thus the mast 
length. For example, if the descending node occurred at 10;00 a.m. instead 
of 9:30 a.m., the array area would become A = 23.78 ra (256 foot ) and the 
required mast length would be increased to 5.33 m (17.5 feet). 

3,5.3 Recommended Solar Array Configuration 

Each of the two solar array configurations — the fixed array and the 
driven array — described above are satisfactory. The fixed array has the 
advantage of not requiring a drive mechanism, but it has the disadvantages 
of increased size, weight and cost. Also, the fixed array requires a much 
longer antenna mast; this extra mast length is undesirable because it 
affects the spacecraft and antenna flexible dynamics. The driven array 
would exhibit dynamic interaction with the spacecraft, but this would 
probably be a smaller effect than that due to the structural flexibility 
of the large fixed array. Based on the above tradeoffs, the driven FRUSA 
array as shown in line 1 of Table 3-14 is recommended. 
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4.0 -DISTURBANCE TORQUES AND MOMENTUM MANAGEMENT FOR LANOSAT-D 


This section summarizes the results of a disturbance torque-momentum 
management investigation for the present concept of the LANDSAT spacecraft 
configuration. The purpose is to' determine magnetic unloading torque 
levels and associated peak stored momentum values which will indicate if 
standard reaction wheels and torquer magnets are sufficient for the design 
of the system. Detailed disturbance torque models and an ideal spacecraft 
with zero attitude errors are used to establish torque and momentum pat- 
terns and magnitudes. Based on the disturbance torque signals for one 
orbits a Fourier coefficient set is compiled which approximates the spec- 
tral content of the disturbances. These may be used in later simulations 
as simplified disturbance models. Eclipse conditions occurring for the 
9:30 a.m. sun-synchronous orbit are included in the disturbance model. 

Based on the results of this study, it appears that the single 50,000 
pole-cm magnet windings are sufficient for momentum unloading. Thus the 
two windings per magnet may be used in a redundant mode. Only a short in- 
terval of magnet command limiting at 50,000 pole-cm is observed per orbit. 

Momentum excursions of less than 2.0 N-m-s occur at all axes. This 
is considerably less than the allowable excursion of 7.32 N-m-s for speed- 
biased operation of the wheels with equal allowed excursions for all 
wheels. Thus it is concluded that the 20 N-m-s NASA Standard Reaction 
Wheels are sufficient for the design. 

Sensitivity to antenna mast heights is not dramatic for the heights 
presently under consideration. These effects, and the effects of antenna 
orientation, are thus not a major concern with regard to momentum storage 
and disturbance torques. 

4.1 Spacecraft Model 

The diagrams in Figure 4-1 display the pertinent measurements for the 
LANDSAT. Since antenna height has not been definitely decided upon, it is 
treated as a parameter to obtain an estimate of its effects on the 
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disturbance torques. Various antenna orientations, primarily looking for- 
ward or zenith, are considered to also determine their effect on distur- 
bance levels. These orientations are chosen because the antenna coverage - 
pattern is approximately a hemisphere over the minus 7 body axis. The 
spacecraft model consists of a standard configuration for the MMS ACS, CDH, 
power modules, and adapter plate. For simplicity, the payload is modelled 
as a cylinder with an end plate since more refined surface models are not 
available. The sun pointing array is modelled as a flat plate that rotates 
at a cone angle of *30 degrees in body coordinates, and the antenna is 
modelled as a thin wide cylinder at mast heights of 3.05 m (10 feet), 

4.57 m (15 feet) and 6.10 m (20 feet). 

Since the array rotates in body coordinates, its center of pressure 
varies for aerodynamic torque calculations. This time varying function is 
modelled in the simulation. The time varying surface normal to. the solar 
panel from the sun is also computed. 

Moments of inertia for the present configurations are I ,, = 2033.7 

on n XX 

kg-m*^, I = 2670 kg-m*^, and I - 2372.7 kg-m^. Weights are 1496.6 kg 

Jj ZZ 

for the spacecraft and 22.7 kg for the solar array and for the antenna. 

The relatively small effects of a time varying inertia matrix are not mod- 
elled. Cross products of inertia are also time-varying, depending primar- 
ily on array position. These affect primarily the gravity gradient torques, 

2 

and are modelled as an average of 271.2 kg-m in all axes. 

For a specified antenna mast size, the center of mass (c.m.) of the 
spacecraft moves slightly as the antenna is rotated. The maximum distance 
is on the order of 1 inch (2.54 cm), and thus a fixed average c.m, is 
chosen for the model. Changes are included for various antenna heights. 

Table 4-1 summarizes the solar-aero surface model data for the present 
study. An antenna height of 4.57 m (180 inches) is shown in the table. 


*Based on Grumman layout drawing; a 9:30 a.m. orbit would correspond to 
an array cant angle of 37.5 degrees. The difference in the resultant 
disturbance torque is insignificant. 
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Table 4-1. LANDSA^Iolar-Aero Model 


Surface Description 

Surface 

Number 

Area 

(1n^) 

Center of Pressure (in) 

Inward Normal Vector 

"x 

m 




"z 

Solar Array Front 

1 


Time Varying, Inernal 

Time Varying, 

Internal 

Solar Array Back 

2 








Interface Back 

3 

828. 

48 

0 

0 

1 

0 

0 

(Spare) 

4 

0 

0 

0 

0 

0 

0 

1 


5 

2304 

24 

31.9 

0 

0 

-1 

0 

ACS Module 

6 


24 

22.9 

-24 

0 

0 

1 


7 

864 

24 

22.9 

24 

0 

0 

-1 


8 

2304 

24 

-15.95 



B 

0.866 

CDH Module 

9 

864 

24 

19.5 



RH 

0.5 


10 

864 

24 

-32.2 

-7.79 


B 

-0.5 


11 

2304 

24 

-15.9 

m 



-0.866 

Power Module 

12 

864 

24 

19.5 

■M 



-0.5 


13 

864 



24 

-32.2 

7.79 



0.5 

Spacecraft Bottom 

B 


0 

0 

0 



0 

Interface Front 

MM 



0 

0 


0 

0 

Payload Front 

■H 



0 

0 


0 

0 

Antenna Front 

17 

B 


-66 


0 

0 

1 

Antenna Back 

18 



-66 


0 

0 

-1 

Payload Body Cyl. 

19 

B 

100 

0 

B 

1 

0 

0 

Antenna Edge 

20 

B 

72.58 

-66 

■■ 

0 

0 

1 


*Dia = 62.2, Length =89.9 
**Dia = 76.0, Length = 10.0 



Acc. Coef. 

0.2 

mm 

0.56 


0.56 


0.56 


B 


■■ 

mgm 

0.56 

0.85 

0.56 

9B 

0.56 

mgm 

0.56 

0.85 

0.56 

0.85 

0.56 

0.85 

0.56 

0.85 

0.56 

mm 

0.56 


0.56 

0.85 

0.56 

BBI 

0.56 

B^B 

0.56 

0.85 

0.56 

0.85 
































































4.2 Disturbance Models 

A standard detailed earth magnetic dipole field model is used as the 
source of the magnetic disturbance torques and for the generation of the 
unloading torques. It is described, for example, in [1] and is not re- 
peated here. A spacecraft residual magnetic moment of 10,000 pole-cm.at a 
fixed random orientation is used. Based on a commonly used estimate of 
1 pole-cm per pound of weight, the 10,000 pole-cm value is conservative. 
Provision is made for rotation of the dipole field as a result of earth 
rotation, but over the short time of one or three orbits considered here, 
this effect is very small and 1s set to zero. 

Solar torques are computed as a function of position in orbit and 
relative sun location using the detailed surface model. An eclipse condi- 
tion is included to zero the solar torques within an angle of + 56.7 de- 
grees of the ascending node*. Shading effects which may be present during 
portions of the orbit are not included. 

Aerodynamic torques are computed with the detailed surface model used 
for the solar torques. Gravity gradient torques are computed using the 
complete theoretical expressions, but the use of an idealized spacecraft 
with zero attitude errors results in constant gravity gradient torques re- 
sulting from the assumed cross products of inertia. Otherwise the peri- 
odic orbital variations in other disturbance torques would result in small 
attitude errors and resulting periodic content in the gravity gradient 
torques. 

The theoretical expression for the gravity gradient torque is given 
by 


Tg = 3 Wq R X (I • R) 


*Assuming an equinox. 
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where 

Wjj - is the orbit rate 

"R - is the distance to the center of the earth 
T - is the inertia dyadic 

From the preceding equation, the gravity gradient torque expressions 
are 

"*^gx ^ ^ '^0 ^^^zz“^yy^ “ ^xy ® ^yz^ 

**^gy * ^ '^0 ^"^^xx-^zz^ ® " ^xy “ ^xz^ 

"*'gz ~ ^ '^0 ^^xz ^ ^yz 

and by setting the pitch angle error e to zero, the following expressions 
are obtained for small (}i angles. 



<!> 


3 w! 


^'^zz"Vy^ 


N-m 

rad 



< 1 ^ 




N-m 

rad 





N-m 

rad 


These may be used to estimate gravity gradient torques for small (]> angles. 
Similarly, the relations 
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N-m 

rad 




<'xx-Izz> 


N-m 

rad 





N-m 

rad 


yield torques for small 9 attitude errors. 

4.3 Momentum Unloading System 

Magnetometers are used to sense the earth's F field, and a magnetic 
moment M is generated by the three-axis arrangement of the bars. The com' 
mand M is gen’erated according to the linear law 


M = -K B X H 
• e 


where K is a suitable constant and is the momentum error vector. The 
preceding unloading law may be arrived at either by minimization of an 
error e between M x F and -K or by geometric considerations. This con- 
trol law is simulated to obtain continuous momentum unloading of the reac- 
tion wheels. 

The dipole magnets for the generation of M are available in standard 
sizes of two 50,000 pole-cm windings on a single core. Those windings may 
be used singly in which case one winding is redundant for backup purposes, 
or they may be used together to obtain 1 x 10 pole-cm with no redundancy. 
The choice is determined by the on-orbit torque levels that are required. 

4.4 Stored Momentum Determination 

It is desired to use the NASA Standard Reaction V/heel (SRW) for reac- 
tion torque and angular momentum storage. This is a 20 N-m-s wheel and in 
the present application four of them will be used in a speed-biased 
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arrangement. It is thus necessary to determine if the orbital momentum 
variations remain below the maximum allowable as measured from the nominal 
speed-biased momentum value. In fact, it is desired to operate the wheels 
as closely as possible to the speed bias. 

To do that, disturbance and magnetic unloading torques are computed 
for an idealized spacecraft in a simulation. The net external torque on 
the spacecraft is transformed to inertial space, integrated to provide in- 
ertial momentum, and transformed back to body coordinates. Since perfect 
attitude control is assumed, the momentum in body coordinates represents 
stored wheel momentum as a function of orbital position. The roll, pitch 
and yaw momenta may be resolved into wheel axes to obtain the excursions 
from nominal speed-biased momenta. 

4.5 Results 

Initial results indicated that magnet signals slightly in excess of 
50,000 pole-cm would be required for momentum unloading. Thus magnet com- 
mand limits of 50,000 pole-cm were applied to determine if unloading could 
be accomplished using just one winding. The results are summarized in 
Table 4-2 for 4.57 and 6.10 meter antenna heights. 


Table 4-2. Maximum Magnet Commands and Reaction Wheel Momenta 


Antenna, Position 


I'’! max N""-" 

4.57 m, looking -2 

61 ,864 

1.78 

4.57 m, looking +X 

60,073 

1.86 

6.10 m, looking +X 

58,947 

1.85 
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The magnet conmands and limiting are more apparent in the plots of 
Figure 4-2a. There it is seen that only the pitch magnet is saturated .at , 
50,000 pole-cm, and only for a brief portion of the orbit near the descend- 
ing node.. Although the desired unloading level is not achieved during 
this time, the system carries the momentum error until it can be unloaded 
thus keeping the resultant momentum less than 2.0 N-m-s. Later paragraphs 
examine allowable momenta levels, but accepting the 1.86 N-m-s as permissi- 
ble it can be concluded that the single winding 50,000 pole-cm is suffi- 
cient. Thus the magnets may be used in redundant mode. 


The total external disturbance torques are shown in Figure 4-2b. 

Since eclipse occurs at + 56.7 degrees from the ascending node, a check is 
easily made on the magnitude of the jump in the solar torques. Without 
allowing for antenna effects, the values should be 


4.07 X 10 
4.07 X 10 


cos 

56.7° 

= 2.23 X 

sin 

56.7° 

= 3.40 X 


10"* 

N-m 

roll 

10-* 

N-m 

yaw 


These agree closely with the plot values. 

From Table 4-2 and from the series of plots in Figure 4-2, it is also 
seen that the sensitivity to antenna mast height and antenna orientation 
is fairly small. This is of interest in estimating effects of possible 
configuration changes. Studies made with a 3 m antenna height show distur- 
bance torques and stored momenta slightly less than for the 4.57 m antenna 
mast. Since the latter are acceptable, the 3 m results are also acceptable 
and are not tabulated here. 


For four reaction wheels in speed-biased operation, with one wheel on 
each of three major body axes and the fourth wheel equiangularly placed 
from the other three, a vector diagram shows that the fourth wheel should 
have time the (equal) momentum in the other three wheels for a net 

zero momentum. And to give all wheels equal momentum ranges, a simple 
calculation shows that the bias momentum should be 7.32 N-m-s for wheels 
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Figure 4-2c, 4.57 m Antenna, zenith Pointing: Wheel Momenta 
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on the axes and -12.68 N-m-s for the fourth wheel. This will allow the 
fourth wheel to accumulate -7.32 N-m-s before momentum saturation, and 
will allow each of the other three wheels to lose 7.32 N-m-s before drop- 
ping to zero speed. 

The preceding bias set point may be tailored slightly to specific or- 
bits if it is known that the stored momentum will show a bias. However, 
the range of + 7.32 N-m-s is used here. Based on it, the 1.86 N-m-s maxi- 
mum excursion shown in Table 4-2 can be easily handled by the proposed 
20 N-m-s SRW's. 

Figure 4-2c contains plots of the roll, pitch, and yaw wheel momenta 
which show that the roll momentum is approximately zero biased. Although 
the initial momentum was zero at the ascending node, a complete orbit shows 
that there is a net accumulated momentum. To insure that the momentum does 
not continue to accumulate on succeeding orbits, three-orbit simulation 
plots are shown in Figure 4-3. These .show more clearly that the roll mo- 
mentum is zero biased and periodic. However, although the pitch and yaw 
momenta are periodic after an initial transient, there is a net bias momen- 
tum. 


This bias is caused by the magnetic unloading law which requires that 
a certain momentum error Hg exists before any unloading torques are gen- 
erated. The observed bias values thus reflect the equilibrium condition 
of the magnetic unloading law. These biases may be reduced by increasing 
the unloading gain K or by changing the magnetic control law to include a 
lead term representing an estimate of the disturbance torque in body co- 
ordinates. In inertial coordinates, the momentum biases are about equal 
in all axes, but in body coordinates they mainly appear in pitch and yaw. 
The results of the present investigation into the adequacy of the wheels 
and the magnetic torque bars are not affected by the biases since they are 
relatively small compared to the reaction wheel momentum capacity. Thus 
reduction of the biases will not be further pursued at this time. 
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Figure 4-3c. 4.57 m Antenna, Pointing Ahead: Wheel Momenta 










For the plots shown in Figures 4-2 and 4-3, the spacecraft was started 
in orbit at the ascending node with zero wheel momenta. Other plots not 
included here show that the wheel momenta appearance varies when the space- 
craft is started at orbit positions of 90 and 180 degrees with zero initial 
momenta. The reason is that an unloading error exists depending on the 
relative orientations of the and B vectors, which in fact depends on the 
initial conditions. Additional studies should examine the possibility of 
reducing the bias momentum on succeeding orbits. This may require a more 
detailed study of unloading gains. 

Figure 4-4 shows the momentum for the three axes as it accumulates 
when no unloading torques are conmanded. The RSS value for one orbit is 
approximately 4.6 N-m-s. In the pitch axis, the aerodynamic torque on the 
antenna (looking ahead) is primarily causing the constant buildup, while 
the small superimposed oscillation is from solar effects. 

To show the sensitivity of the disturbance torques to the mast height 
and antenna orientation, the Fourier coefficients for the fundamental and 
first harmonic are tabulated in Table 4-3. These may also be used later 
for simulation models. In the Table, T , T , T denote the solar roll, 

^ O Jr o ^ 

pitch, and yaw torques, T,^, T,, . T^,, the aerodynamic roll, pitch, and 

aX fiy aZ 

yaw torques, T , T , T , the gravity gradient roll, pitch, and yaw 
Qa gy gz 

torques, Tj^^, T^^^, the residual magnetic torques, and T^^, T^^, T^^ 
denote the total disturbance torques. The Fourier model is 

2 

‘''ij " " "o* 

n=l 


where 


i = s, a, g, m, d 
a = X, y, z 
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Table 4-3. Fourier Coefficients for Disturbance Torques 


4.57 m Antenna, Pointing Zenith 


Torques 

Fourier Coefficients 






N-tn 




^2 

4 

Tsx 

-12.5 X 10"^ 

2.19 X 10"^ 

-2.17 X 10"^ 

-1.03 X 10"^ 

1.20 X 10'^ 


4.07 X 10“^ 

-6.49 X 10“^ 

-1.05 X 10"^ 

2.39 X 10"^ 

3.09 X 10’^ 


-4.66 X 10'^ 

2.63 X 10"^ 

-2.97 X 10"^ 

3.32 X 10’^ 

1.09 X 10’^ 

Tax 

2.77 X 10“^ 

-1.26 X 10"^ 

-1.15 X 10'^ 

1.00 X 10'^ 

6.07 X 10’^ 

Tay 

-2.96 X 10“^ 

1.53 X 10“^ 

2.36 X 10’^ 

3.88 X 10'^ 

-2.49 X 10“^ 


-4.68 X 

-3.29 X 10"^ 

-3.00 X 10"^ 

3.15 X 10‘^ 

-6.72 X 10”^ 

Tgx 

9.10 X 10"^ • 

0. 

0. 

0. 

0. 

V 

-9.10 X 10'^ 

0. 

0. 

0. 

0. 

Tga 

-1.53 X 10'^^ 

0. 

0. 

0.- 

0. 

*^mx 

0. 

0. 

0. 

0. 

0. 


-2.10 X 

2.64 X 10"® 

-4.47 X 10‘^° 

-4.19 X 10'^° 

-1.14 X 10“^ 


1.87 X 10"^ 

0. 

0. 

0. 

0. 

^dx 

7.95 X 10"^ 

2.18 X 10“^ 

-3.32 X 10'^ 

-9.34 X 10“^ 

7.27 X 10"^ 

^dy 

-8.98 X 10"^ 

8.79 X 10‘^ 

-4.34 X 10“^ 

6.26 X 10'^ 

-2.18 X 10"^ 

Tdz 

-4.81 X 10‘^ 

2.30 X 10"^ 

-3.27 X 10’^. 

3.47 X 10“^ 

4.15 X 10'^ 
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Table 4-3. Fourier Coefficients for Disturbance Torques (continued) 

4.57 m Antenna, Pointing Ahead 


Torques 


Fourier Coefficients 


N-m 

a 

0 

"l 


^2 

■^2 


-1.19 X 10“^ 

2.10 X 10"^ 

-2.17 X 10'^ 

-1.04 X 10"^ 

1.28 X 10'^ 

sx 




. 



3.97 X 10'^ 

-6.22 X 10"^ 

2.09 X 10'^ 

2.35 X 10'^ 

-3.24 X 10“^ 


-5.07 X 10“.^ 

2.75 X 10'^ 

-3.02 X 10"^ 

3.84 X 10'^ 

1.08 X 10"^ 


2.77 X 10'^ 

-1.26 X 10"^ 

-1.15 X 10'^ 

1.00 X 10"^ 

6.07 X 10~^ 

^ay 

1.49 X lO""^ 

1.53 X 10'^ 

-5 

2.36 X 10 

3.88 X 10"^ 

-2.49 X 10'^ 


-5.42 X 10"^ 

-3.29 X 10"^ 

-3.00 X 10”^ 

3.15 X 10"^ 

-6.72 X 10"^ 

T 

9.10 X 10"^ 

0. 

0. 

0. 

0. 

gx 






V 

-9.10 X 10"^ 

0. 

0. 

0. 

0. 


-1.53 x00“^^ 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

mx 






T 

my 

2.10 X 10'^° 

2.64 X 10"^ 

-4.47 X 10'^° 

-4.19 X 10'^° 

-1.14 X 10"® 

■^mz 

3.23 X 10’^ 

0. 

0. 

0. 

0. 

^dx - 

7.93 X 10“^ 

2.13 X 10"^ 

-3.32 X 10'^ 

-9,40 X 10"^ 

7.35 X 10'^ 

^dy 

-7.20 X 10'^ 

9.06 X lO”^ 

-4.03 X 10"^ 

6.22 X 10"^ 

-2.82 X 10”^ 


-5.60 X 10“^ 

2.43 X 10"^ 

-3.32 X 10"'^ 

3.53 X 10"^ 

4.11 X 10'^ 
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Table 4-3. Fourier Coefficients for Disturbance Torques (continued) 


6.10 tn Antenna, Pointing Ahead 


1 

1 Torques 

N-m 

Fourier Coefficients 


®o 

"l 

*>1 ■ 

- ^2 

^2 

^sx 

"sy 

Tax 

"*^ay 

"az 

^gx 

\x 

T 

my 

\z 

-T 

'dx 

^dy 

-1.17 x 10'^ 

3.95 X 10"^ 
-5.14 X 10'^ 

2.77 X 10'^ 

2.17 X 10"^ 

-5.42 X lO'^ 

9.10 X 10'^ 

-9.10 X 10'^ 
-13 

-1.53 X 10 
0. 

-2.10 X 10"^^ 

3.23 X 10"^ 

7.96 X lO"^ 
-6.54 X 10"^ 
-5.60 X lO"^ 

2.12 X 10"^ 
-6,13 X 10"^ 
2.75 X 10’^ 
-1.30 X 10“^ 
1.53 X 10"^ 
-3,29 X 10“® 
0. 

0. 

0. 

0. 

2.64.x 10"^ 

0. 

2.12 X 10"^ 
9.17 X 10"^ 
2.43 X 10'^ 

-2.17 X 10'^ 
3.40 X 10“^ 
-3.02 X 10"^ 
-1.18 X 10"^ 
2.36 X 10'^ 
-3.00 X 10"® 
0. 

0. 

0. 

0. 

-4.47 X 10"^ 

0. 

-3.36 X 10"® 
-3.90 X 10“^ 
-3.32 X 10"'^ 

-1.06 X 10"^ 
2.30 X 10"® 
3.84 X 10"® 

1.00 X 10"® 

4.01 X 10"® 
3.15 X 10"^ 
0. 

0. 

0. - 
0. 

-4.19 X 10"^° 

0. 

-9.59 X 10"® 
6.32 X 10"® 
3.53 X 10"^ 

1.30 X 10~® 
-6.87 X 10"® 
1.08 X 10"^ 
6.07 X 10“® 
-2.52 X 10"® 
-6.72 X 10“® 
0. 

0. 

0. 

0, 

-1,14 X 10"®’ 
0. 

7.38 X 10"® 
-3.20 X 10"® 
4.11 X 10"® 
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Note that the cross products of inertia contribute small constant 

gravity gradient torques. These may be computed to be 9.1 x 10~^ and 
“4 

-9.1 X 10 N-m in roll and pitch, respectively. These compare very well 
with the table values for roll and pitch. 
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5.0 INITIALIZATION OF FINE CONTROL 


After completion of the initial orbit acquisition, and later on if 
necessary after stationkeeping or any other maneuvers, the fine control 
mode must be initiated. The TORS antenna must be deployed and the solar 
array drive placed into normal operation. 

Normal on-orbit attitude control of the L-D is accomplished using 
reaction wheels (RW's) and only two areas important to fine control mode 
initialization must be examined. 

(1) When entering the RW attitude control mode, the spacecraft body 
rates must be compatible with the RW capabilities and the mini- 
mum impulse bit of the reaction control system (RCS) used for 
momentum unloading. 

(2) An accurate and stable attitude reference for earth pointing 
must be available. 

Since the L-D mission will be flown on an MMS bus, point 1 above hardly 
needs any further examination. The mass expulsion RCS and the torque and 
momentum capabilities of the RW's specified for the MMS [1,2], are com- 
patible, When under mass expulsion attitude control, the MMS is capable 
of reducing spacecraft rates to less than 0.1 deg/sec per axis. The RW 
attitude control system (ACS) is capable of stabilizing the spacecraft 
within 30 minutes from initial rates of up to 0.25 deg/sec about any 
axis [2]. The minimum torque impulse bit available from the 0.2 Ibf 
thruster is conservatively estimated for L-D as 0.05 ft-lb-sec which is 
more than two orders of magnitude smaller than the wheel capacity of 
15 ft-lb-sec {20 Nms) and the ability to unload the wheel is assured. In 
conclusion, there exists no problem in transferring from RCS to RW atti- 
tude control. 

Establishing an accurate and stable attitude reference for earth 
pointing Is thus the main problem encountered when initializing the fine 
control mode and this section is therefore mainly concerned with the 
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nitialization of the MMS stellar-inertial attitude reference system. The 
atter task separates into two sub-tasks. Following the normal HMS acqui- 
iition procedure (described below) it must first be shown that stars 
lighted' by the on-board star trackers can be identified, and secondly, 

:hat the on-board attitude reference algorithm (filter) will converge from 
;he attendant, relatively large, initial attitude uncertainty. Sec- 
tions 5.2 and 5.3, respectively, will address two different methods of 
>tar identification for initializing the on-board stellar-inertial atti- 
tude reference system (ARS). The method presented in Section 5.2 utilizes 
sarth magnetic field measurements to aid in the initial star identifica- 
tion, while the method in Section 5.3 correlates observed star sightings 
vith entries in a reduced star catalog. The convergence of the on-board 
filtering algorithm is demonstrated in Section 5.4. Initialization of 
the fine control mode appears to pose no problems for L-D. 

Once a precise inertial attitude reference of the spacecraft has been 
established, L-D ephemeris data is used to convert the inertial reference 
to an earth pointing reference. Then ACS commands are issued that cause 
the spacecraft to maneuver from an inertial hold mode to its normal 
on-orbit attitude: yaw (Zj^) to nadir, pitch (y^^) normal to the orbit 

plane, and roll (x^) in direction of orbital motion. The RW's will be 
used to perform this maneuver. Any further maneuvers performed during the 
life of the spacecraft, such as stationkeeping, for example, will not 
require a reinitialization of the ARS. At the end of even the longest of 
any such maneuvers, the inertial reference unit (gyros) provides an atti- 
tude reference that will certainly be accurate enough to identify the 
stars sighted fay the trackers. Attitude and gyro bias updates can then 
be computed. 

The TORS antenna should be deployed once coarse sun acquisition has 
been completed. Since MMS/L-D will probably acquire the sun with the 

A 

negative yaw (-2j^) axis, the antenna dish could find itself looking 
directly into the sun, a situation rarely encountered during normal 
on-orbit operations. To assure thermal loads compatible with normal 
on-orbit conditions the antenna elevation angle- should be conmanded to its 


5-2 



maximum dip angle of -26° and the azimuth to +90° so that the antenna 
looks away from the sun and in the positive pitch direction. 

The solar array will be placed into normal operation (see Section 3.5) 
once the spacecraft has maneuvered to its nominal earth pointing orienta- 
tion. Correct initial positioning of the array will be accomplished by 
ground command. Spacecraft ephemeris in ECI coordinates and time of year 
are sufficient information to compute the correct array angle relative to 
the solar array/spacecraft scribe-mark. The correct array orientation can 
be verified from telemetry data of the coarse sun sensor which is mounted 
on the solar array. 

5.1 Nonnal MMS/L-D Acquisition Sequence 

During launch the ACS is in a standby condition and the actuators are 
disabled. After separation from the launch vehicle {shuttle or Delta 3910 
booster) the ACS is enabled and the solar array is deployed and subse- 
quently rotated to a predetermined position relative to its scribe-mark. 

For L-D the solar cells will look toward -z^^. The nonnal MMS acquisition 
sequence is to be autonomous by utilizing the on-board computer. 

Coarse Sun Acquisition 

After array deployment the spacecraft acquires the sun with its nega- 
tive yaw axis {-z^ axis) and the spacecraft rates are reduced to an appro- 
priate level. This maneuver is accomplished by the on-board computer using 
coarse sun sensors located on the solar array, the Inertial Reference Unit 
(gyros) for rate information, and the reaction control jets, magnetic tor- 
quers and/or reaction wheels to provide control torques. Normally coarse 
sun acquisition seeks an orientation with the solar array normal to the 
sun line. Because of the array cant angle of 37.5° (9r30 am orbit) the 
sun line would then not be aligned with one of the principal control axes 
of the spacecraft. Since it may be required to perform a rotation about 
the sun line to acquire appropriate stars into the fields of view of the ‘ 

A 

star trackers, it is preferred to align the ~Z|^ axis with the sun which 
keeps the solar array at 52.5°(90°-37.5°) to the sun line. The 20 percent 
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loss in generated electrical power is not critical since the payload 
sensors, which consume more than 50 percent of the power, are not yet 
activated. For one axis^the output signals of the array-mounted coarse 
sun sensors must then either be biased by 37,5“ to account for the array 
cant angle, or the coarse sun sensors must be mounted to the uncanted 
portion of the solar array shaft. 

Utilizing the mass expulsion system, the ACS is capable of completing 
a coarse sun acquisition sequence within 15 minutes in a sunlit segment of 
the orbit for initial rates of up to 2 deg/sec about each control axis. 

This would apply to acquisition after separation from a Delta 3910 launch 
vehicle. Utilizing reaction wheels, which would apply to acquisition 
after a shuttle launch, the ACS can complete coarse sun acquisition within 
30 minutes in a sunlit segment of the orbit for initial rates of up to 
0.25 deg/sec about each control axis. Coarse sun acquisition is con- 
sidered complete when the attitude error between the sun line and the 
coarse sun sensor reference axis (biased) is less than ten (10) degrees 
(including albedo effects), the maximum rate about each of the control 
axes is less than 0.05 degree per second and control is maintained by the 
reaction wheels. Following coarse sun acquisition, the magnetic torquers 
and/or the mass expulsion system are used to reduce the reaction wheel speeds 
to less than 25 percent of no load speed. 

Fine Sun Acquisition Node 

In the fine sun acquisition mode the coarse sun sensor is replaced by 
the fine digital sun sensor located in the ACS module with reference axis 
aligned with the axis. The initial conditions for the fine sun acqui- 
sition mode are defined by the final conditions of the coarse sun acqui- 
sition mode. Fine sun acquisition shall be considered complete when the 
attitude error between the sun line and the fine sun sensor reference is 
less then 0.1 degree, the maximum rate about each of the control axes is 
less than 0.05 degree per second, and the reaction wheel speeds are less 
than 25 percent of no load speed. Fine sun acquisition can be completed 
vithin ten (10) minutes of entering the mode. 
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After completion of the fine sun acquisition mode the spacecraft 
enters the stellar acquisition mode, which, as was mentioned above, is 
part of the main topic of this section. 

Stellar Acquisition Mode 

The initial conditions for the stellar acquisition mode are defined 
by the fine sun .acquisition mode. Stellar acquisition shall be considered 
complete when the attitude of the ACS reference (optical cube) in the ACS 
module is aligned within 0,1 degree (each of three axes) of a specified 
inertial attitude (including verification of required guide stars), the 
maximum rate about each of the control axes is less than 0.01 degree per 
second, and the reaction wheel speeds are, less than 25 percent of no load 
speed. Stellar acquisition must be completed within four hours after 
entering the mode. At the end of stellar acquisition for spacecraft is 
considered to be in an inertial hold mode from which its normal on-orbit 
orientation is acquired by performing the appropriate maneuver. 

5,2 Star Identification/Attitude Determination From Magnetic Field 

Measurements 

Conditions at the completion of fine sun acquisition establish the 
initial conditions for the start of attitude determination from magnetom- 
eter measurements. The conditions are that the negative yaw axis is 
within 0.1 degree of the sun line and that the maximum rates are 
0.05 deg/sec about any of the control axes. The attitude about the sun 
line is unknown and is to be estimated. The preferred way to do this for 
low orbit missions is to wait for a point in the sunlit portion of the 
orbit at which the earth's magnetic field is oriented more or less ortho- 
gonal to the sunline (60° or more). For high inclination missions, such 
as L-D, the portion of the orbit around the poles is generally not satis- 
factory, One purpose here is to point out angular error sources in the 
overall magnetic attitude determination scheme. 

5.2.1 Magnetic "B Field Modeling 

Autonomous attitude determination is desired and to do this a suffi- 
ciently accurate magnetic B field model must be available on-board at the 
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time of attitude determination. Comparisons of a simple tilted dipole 
field model and the standard spherical harmonic model reveal that the 
"B field vectors computed by the models may differ in direction by as Tnuch 
as 10-15° or more. This rules out the possibility of using the tilted 
dipole model on-board for autonomous attitude determination since the star 
tracker field-of-view is only 8° x 8° and a 10° - 15° uncertainty would 
also be unsatisfactory for star identification. Since launch date and 
therefore approximate vehicle location in the B field is known before 
launch, a local portion of the spherical harmonic model may be stored into 
the on-board processor. This would still permit autonomous attitude 
determination. However, storage of a complete harmonic model may be pro- 
hibitive, and if it should be required, the model would have to be com- 
puted on the ground and sent up for use at the proper time. 

Finally, the agreement between the harmonic model and nature has to 
be determined since it is a direct error source in the attitude determina- 
tion. Estimates of the divergence between the hannonic model and nature 
are not well known because of a lack of data, but published results [7] 
show that solar-magnetic field attitude determination produces accuracies 
of 1 to 3.5 degrees from field measurements that are in error by up to 
9.2 degrees. Repeated measurements at successive points in orbit do, how- 
ever, help to reduce the field uncertainty effects. 


To be of use in the attitude estimation process, the magnetic field 
must be available in a proper coordinate frame for on-board use. This 
coordinate frame is referred to as the sun frame and is defined in detail 
in the next section. Typically, the tenth order hannonic model B field 
is computed in spherical coordinates as 


B = 


B 

B 

B 


r 

e 


where r is the radius from the geocenter to the field point, A is the East 
longitude, and e is the colatitude. This field vector is then transformed 
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suitable routine orthogonal transformations to obtain it finally in 
sun coordinates. Possible errors in the final field vector as a result 
of unknown angular errors in the transformations are discussed after the 
coordinate systems are defined. In any event, the magnetic field is 
assumed to be transformed error- free into a standard set of Earth-Centered- 
Inertial coordinates which will be related to the sun coordinates by 
standard orbit transformations. The error-free assumption is justified 
because up to that point the transformation does not involve orbital 
geometry. 

5.2.2 Coordinate Transformations and Frames 

It is convenient to have a set of Earth-Centered-Inertial coordinates 
Xj yj Zj with Xj toward the vernal equinox point, Zj normal to the 
equatorial plane, and y^ completing the set. It is in this set that the 
magnetic field is assumed available by error-free transformation from 
magnetic coordinates. 

To perform the attitude determination, it is convenient to establish 
a set of sun coordinates x^, y^, centered in the spacecraft with 
-Zg toward the sun, in tne orbital plane, and y^ completing the set. 

Then the vehicle body coordinate set Xg, yg, Zg is defined by the respec- 
tive rotations a, g, y where a is about the x^ axis, g is about the inter- 
mediate y axis, and y is about the final Zg axis. Angles a and g are less 
than or equal to 0.1 degree as noted earlier, and y is arbitrary about the 
Zg axis which is nominally the anti -sun line. Figure 5-1 shows the frames 
and angles involved. 




The transformation from sun to body coordinates is [A]gg given by 

cosy cosg COSY sing sina + siny cosa -cosy sing cosa + siny sina 

✓ 

-siny cosg -siny sing sina + cosy cosa siny sing cosa + cosy sina 
sing -cosg sina cosg cosa 


and is explicitly required.for the least-squares attitude determination in 
Section 5.2.3. 
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Additional transformations are required also, but those involve angu- 
lar quantities such as orbit inclination, orbit position angle, and angle 
of the ascending node which are all known quite accurately. For instance, 
spacecraft position is assumed known to within 100 meters which, at an 
altitude of 705 km (388 nautical miles) creates an angular error of 
approximately 3 arc-sec. Errors such as this in the computation of trans- 
formations are entirely negligible compared to the potential disagreement 
between the magnetic model and nature so they may be ignored. The errors 
are also negigibly small compared to the 8° x 8® field'-of-view for the 
star tracker. Thus the overall transformation from ECI coordinates to sun 
coordinates is designated as and with it the magnetic field avail- 

able in ECI coordinates may be transformed to the sun coordinates by 

Bs = Wsi^I 

for use in the attitude determination. 

5.2.3 Attitude Determination 


Magnetic field measurements are obtained from the magnetometers in 
body axes, and since the field model is known in sun coordinates as 
the measurement equation may be written as 


= [A] 


BS 


where z is the vector of three measurements 


xs 


ys 


B 


zs 


The additive noise is represented by the vector V, and since the measure- 
ments are processed at discrete times in orbit, the covariance matrix of 
V may be written as 
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Erv(t,l = [R] 6.J 

where [R] is a constant matrix. Modeling the noise as uncorrelated between 
the axes and with equal variance in all axes permits writing [R] as 


[R] 


r 2 

% 

0 

0 


0 

2 

% 

0 


0 

0 


For the magnetometers considered here, the resolution is 0.004 gauss 
which yields a = 0.001 gauss if a rounding-error model with a uniform 
distribution is used. The noise is considered zero mean although bias 
may result since the accuracy of the magnetometers is +1% of the full 
scale value of jt).5 gauss. 


The emphasis here is on a computationally simple algorithm that will 
produce estimates for the angle about the sun line accurately enough so 
that unambigous star identifications can be made and the stellar-inertial 
attitude reference system can be initialized. Thus the extended Kalman 
filter is not considered for magnetic attitude determination. Also of 
importance is the fact that rates about the sun line can be appreciable 
at this time and thus the yaw attitude may change betv/een successive 
magnetometer attitude determinations. Thus the yaw gyro rate has to be 
integrated so that the yaw angle change can be propagated between succes- 
sive attitude determinations. Since the yaw axis is within 0.1 degree of 
the sunline, a and e are always small and accurate enough for stellar con- 
vergence; however they may also readily be estimated by the least squares 
procedure if desired. The procedure is as follows. 

Let x° represent the initial best estimate of the vector x = 

Since the measurement equation is nonlinear and implicitly involves x, 
it may be written as 


z = h(x) + V 
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h(x) 


• cnc6B^^ 

+ (cySBSa + SYCa) 

+ 

(-CysBca + sysa) B 

-sycsB^j 

+ (-SySBSa + CyCa) B 

ys 

■+ 

(sySBca + cysa) 


- cBsa B 

+ 

cecc. B^j 


and linearized about x = 


{0 , 0 , Y°) to obtain the matrix H defined as 


H 



siny° B ^3 

-COSy° 

cosy° B ^3 

siny° 

- B 

B 

ys 

xs 


-siny 

-cosy - Siny 


The complete linearization may now be written as 


AZ = H Ax + V 

where the x and z are defined as 


- ~ **o 

AX = X - X 

- 2 - 


These equations would be suitable for estimating the attitude x based on 
Ax computed from the relation 

AX = az 


However, in the present case a and 6 are known within 0.1° if the fine sun 
sensor is operative. This accuracy is probably greater than can be 
obtained from the magnetic field, and thus a and B do not normally have 
to be estimated. The preceding formulation is still applicable with some 
notational changes as follows. 
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Using only the last column of matrix H, define 


M = 


- sir>Y 

- - sinv° 


Then the measurement equation may be condensed to 


AZ^ 


■''r 


= M Ay + 


CVI 

ISI 

< 

1 




from which a least squares estimate for Ay may be obtained as 

r. 


Ay = (h"'’ M)"^ m"*" 


AZ 


1 


Mr 


The attitude yis then obtained from- 

Y - V° + (iJ - [A'(?)g 5 ] i^) 


( 5 . 2 - 1 ) 


where Az is understood to be 


Az = 


AZn 


AZ, 


and the primed transformation A'(x°)gg is obtained from A(x®)g2 by delet- 
ing the last row, as 


CA'] 


BS 


cosycose cosysinBsina + sinycosa 
-sinycosg -sinysinBsina + cosycosa 


- cosysinecosa+ sinysina 
sinysinBcosa + cosysincc 
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The first estimate of 7 ° will be determined shortly, and then B , 

[A'(x )]gg, and the measurement i(t^} are all known so that y may be com- 
puted. Let the first estimate be y(t-j), and determine the incremental 
attitude change y^ up to time b^ by integration of the gyro outputs. Then 
the new initial estimate for the second least-squares determination is 
given by 


y°(t2) = y(t^)+yg(t2) 

Using this y°(t 2 ), the computation in (5.2-1) may be repeated to determine 
the second attitude estimate y(t 2 ). The process is repeated until y(t^) 
has converged, i.e., successive computed updates Ay are small and have 
zero mean. 

The first estimate x°(t^), i.e., y'^(t-j), may be determined as follows. 
Since a and p are within 0.1 degree, they may be set to zero in the 
measurement equation 


z = [A]gg 


Writing the equation out and solving for the remaining unknown y° yields 
the expression 


= cos 


-1 


+ Zo B 


-1 ^sx 


2 ^sy 




sy 


which provides the starting estimate. 

The flow diagram v^hich follows contains the main steps for the 
magnetometer attitude estimation. 
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Process first B measurement to get 
x°(t,) =■ (0.0, 

and [A’(i“(t,))]gj 


Initialize integration of gyro 
outputs to begin obtaining 
incremental attitude change Yn(to) 


Evaluate the attitude estimate 
Y(t^) = T°(t^) + M‘^(z(t^) - AB^) 


Propagate the attitude by 
integration of gyro outputs 
to obtain Yg(^2^ 


From the new initial estimate 

Y (^2 ® Y{i^]) '*■ Yg(t2)' 

Take the magnetometer reading 
at time to 


Compute the new 


Reinitialize integration 
of gyro outputs 


Evaluate the attitude estimate 
Y(t2).= Y°(t2) + (M’^M)'^M’^(i(t2) - A'B^ 


Return and 
Increment 
the time step 


Begin Stellar Attitude 
Acquisition Procedure 


Figure 5-la. Attitude Processing Flow Chart 
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The arrangement of the preceding flow chart is arbitrary in several 
places, and it is implicitly assumed that the computational time per esti- 
mate is negligible compared to the time difference t. - t. between esti- 
mates. 

5.2.4 Error Sources and Effects 

This section summarizes the potential error sources and indicates the 
accuracy expected for the magnetic attitude determination. The error 
sources and their effects are: 

(1) Divergence between the magnetic model and nature. A mean 
angular error is 4.3° extracted from Reference 1, with a 
maximum error of 9.2°. Effects of this divergence are reduced 
by repeated measurements at several points along the orbit. 

(2) Magnetometer measurement noise of 0.004 gauss. These effects 
are also reduced by repeated measurements along the orbit. 

(3) Errors in coordinate transformations. Since the transforma- 
tions involve angular quantities that are accurately known, 
the error in attitude determination that results from an 
inaccurately transformed B field is negligible. 

Results in Reference 7 show that attitude errors on the order of 
1.56° (average) with a maximum of 3.5° are achievable. Such errors are 
well within the 8° x 8° field-of-view of the star tracker, and thus the 
solar-magnetic determination should be accurate enough to permit initiation 
of stellar acquisition. 

5. 3 Star Identification By Correlation 

Another method of determining the orientation of the spacecraft about 
the sun line is by correlating star tracker observations with a star 
catalog. The star catalog contains only those stars above a certain 
brightness which should be encountered by the spacecraft's star tracker 
as it sweeps out at a swath about the sun line. The correlation method 
will now be briefly described. 
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5.3.1 Description of Method 

-The spacecraft Is commanded to rotate about the sun line at a fixed 
rate, say 0.1 deg/sec, and the star tracker sensitivity threshold is set 
such that it will detect only stars above a specified brightness which 
matches the brightness of the stars included in the catalog. At the 
beginning of the spacecraft rotation, or scan, the inertial reference sys- 
tem is initialized with a guess of zero rotation angle about the sun 
line, the rotation angle being referenced to a suitably defined sun point- 
ing reference frame. The stars in the star catalog are sequentially 
ordered according to this reference point and each star is specified by an 
azimuth or scan angle measured from the reference (in the direction of 
spacecraft rotation) and a declination relative to the swath plane, <5^ 

(the subscript c stands for catalog). ' Then, as the star tracker sweeps 
out the swath, the time and the scan angle for each star encountered by the 
tracker will be recorded. The scan angle is obtained from the inertial 
reference unit, i.e,, integrated gyro rate, and is denoted by <j)j (sub- 
script T for tracker). If the spacecraft is controlled so that it remains 
well aligned with the sunline, one can start to correlate the sequence of 
star observations to the catalog entries after one revolution. If the 
spacecraft rotation axis does not remain well aligned with the sunline, 
some stars may be missed due to the resulting wobble, and 2-3 revolutions 
of data taking may be required. LANDSAT D will be maintained sun point- 
ing to within 0.1 degree when using the fine digital sun sensor or an 
attitude reference and one revolution should suffice. 

The correlation proceeds then as follows. The first star encountered 
by the tracker is assumed to be the first star in the catalog and the dif- 
ference of the scan angles is computed, 

4>q(1) = “ ^T^^^ 

The angle <I>q( 1) may be considered the first update on the initial zero 
estimate for the orientation angle of the spacecraft about the sun line 
and it is added to all the observed scan angles (J)j(n). All the observed 
stars are now correlated with the entries in the star catalog, considering 


5-16 



all possible combinations, and the number of successful correlations or 
matches are recorded. A successful match is considered to occur between 
two stars when the observed star scan angle li)j(n) falls within a specified 
e of the catalog scan angle 4i^(m), i.e. , 

- <{)-j-(n)j < e 

The process is then repeated for the next pair of stars, for example the 
first star encountered by the tracker is assumed to be star No. 2 in the 
catalog and the number of successful correlations are determined under 
this assumption. This process is repeated until all possible star pairs 
have been exhausted, and the pairing that produced the most matches is 
declared the winner and the corresponding scan phasing angle is the 
final estimate of the orientation angle of the spacecraft about the sun 
line. ‘}>Q is not selected as the angle computed for the first forced pair 
in the sequence, but once good correlation for a particular sequence has 
been established, all star measurements are considered and a least squares 
fit to can be obtained. Initial calibration of the scan axis gyro bias 
can also be obtained. In principle, the elapsed time between two star 
observations, the Acjiy scan angle as obtained from the gyros, and the 
for the two stars as obtained from the catalog, contain all the informa- 
tion to compute gyro bias. Namely 


Af _ - A<t>-p 

*>9 = 

^scan 

In practice one uses all the' data and computes a least squares fit estimate 
for the gyro bias. 

Using this method the spacecraft orientation about the sun line can 
usually be determined within about 0,2 degree and the gyro bias estimate 
within 0.2 deg/hr. The errors are mainly due to bias and noise imperfec- 
tions in the tracker and gyros. Residual gyro bias appears to be the 
main error source in the estimate of the scan phasing angle (J)^. The num- 
ber of stars also influences the accuracy; the more data, the better the 
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estimate. The latter is, however, of little concern to LANDSAT D since 
one is not so much interested in the accuracy during initialization, -as 
one is Interested in a fast correlation process using only a few good 
stars. 

A potential disadvantage of the correlation method is the length- of 
time that may be required to obtain enough stars since trackers will not 
read accurately above a certain angular rate which limits the speed of 
commanded rotation about the sun line. Depending on the number of stars 
in the swath/catalog, a more serious requirement is the rather large 
amount of computation required to perform the correlation tests. Since 
the correlation is nominally only performed at the beginning of the mission 
the software, which includes the swath catalog, should not be left in the 
computer. But erasing it and reading a new software from a special 
memory would require a ground command which does not fit the desired 
autonamous operation plan. Also, the ground command may have to be 
delayed until a certain favorable orbit position is reached. Since 
autonomous operation seems hard to attain with this method, it is recom- 
mended to perform the correlation task on the ground. In this case the 
star tracker data is telemetered to the ground and after convergence of 
the correlation procedure, the correct scan phasing angle about the sun 
line and the initial gyro bias calibration results are sent back to the 
spacecraft. 

5.3.2 Simulated Star Identification by Correlation 

Star identification by correlation has been simulated for L-D by a 
digital program using real star locations as obtained from the Yale star 
catalog. The time of year was assumed to be around July 20th such that 
the right ascension of the sun was exactly 120 degrees. The negative yaw 
axis of L-D was pointed to the sun with an accuracy of 0.1 degree and the 
spacecraft rotated about the sun line at 0.1 deg/sec. The cross-axis 
rates were assumed in a sinusoidal limit cycle with peak rates of 
0.05 deg/sec and position amplitude limited to 0.1 degree. Initial uncali- 
brated gyro biases were set at values averaging to 1.5 deg/hr and the 
sensitivity threshold of the S'* x 8° FOV star tracker (BBRC CT-401) was set 
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to My = +4. The latter will result in a small star catalog for the swath 
and provide for faster and less ambiguous star identification/correlation 
thereby taking full advantage of the relatively large field-of-view (FOV) 
of the tracker. 

The tracker slant angle <}>g relative to the ACS module, see Figure 2-5, 
has been assumed as = 30”, and the angle between the lines of sight of 

the two trackers as 70°, i.e., + 35° from the projection of the pitch axis 
onto the plane defined by the lines of sight of the trackers. Since both 
trackers make the same angle of 65.8° with the negative yaw axis, only 
one tracker is used. 

Figure 5-2 shows a plot of all stars of 4th order magnitude and 
brighter and also shows the line traced out by the tracker boresight on 
the celestial sphere. Figure 5-3 removes all stars which do not fall into 
the swath traced out by the field-of-view of the tracker and marks those 
which are eclipsed by the earth. As can be seen, only four stars are 
left making correlation of tracker observations with the catalog very 
simple. Figure 5-3 also shows the assumed starting point of the scan 
according to which the sequence of stars in the star catalog is organized, 
and also shows the actual starting point. Thus, the first star encountered 
by the tracker is actually the star with Yale number 3487 but the tracker 
assumed at first that it was star number 1038. During the correlation 
process all stars are eventually correctly identified and a scan phasing 
angle of = 224.9 degrees is determined. This is within 0.1 degree of 
the actual scan phasing angle of (fi^ = 225 degrees between the initially 
assumed scan starting point and the actual starting point. The yaw gyro 
bias was estimated as 1.2 deg/hr and its actual value was 1.4 deg/hr. 

It is seen that the accuracy of the correlation method is much better 
than determination if the orientation about the sun line using magnetic 
field measurements. However, as shall be shown in the following section, 
this improved accuracy is not required to initialize the stellar-inertial 
attitude reference system and the magnetic method will work just fine. 
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5.4 Convergence of On-Board Attitude Reference Filter 

Attitude of the spacecraft about the sun line was estimated to within . 
about 2° when magnetic field measurements were used, and to within about 
0.1® when the star catalog correlation method was used. The magnetic 
method is simpler and faster and is preferred. 

To demonstrate convergence of the on-board attitude reference filter 
with initial conditions corresponding to the nominal conditions at the end 
of magnetic-solar acquisition, an attitude estimation simulation was per- 
formed using an extended Kalman filter. For the spacecraft this is a 
part of the stellar acquisition mode. Since the negative yaw axis is sun 
pointing to within 0.1 degree, roll and pitch attitude uncertainties are 
a maximum of 0.1 degree at the initiation of fine stellar acquisition. 

These errors along with a nominal 2.0 degree estimation error about yaw 
resulting from the magnetic-solar attitude determination are used as 
initial conditions for the filter. 

Actual stars that the star tracker would see are determined by select- 
ing a time of year for the sun position. As in Section 5.3, the time 
around the 20th of July is again selected with the sun angle at 120° 
measured in the ecliptic plane. As described earlier, when the spacecraft 
is sun pointing the star trackers are inclined at an angle of 65.82° to 
the sun line. The star tracker sensitivity threshold has been set to 
* 5 so that the tracker will sweep out a swath of fifth magnitude or 
brighter stars as shown in Figure 5.4 when the spacecraft is commanded to 
rotate about the sun line. The field-of-view is 8° x 8°, and stars 
eclipsed by the earth are not shown. Also several other stars satisfying 
the My = 5 or brighter criterion, are not shown and the stars indicated 
in Figure 5-4 represent the reference star catalog to be used in stellar 
convergence. The stars are listed in Table 5-1. 

With a 2° attitude uncertainty about yaw at the end of solar-magnetic 
acquisition, it is doubtful whether the stars that happen to appear in the 
field-of-view of the trackers at that time are sufficiently separated from 
neighboring stars to yield a^ unambiguous star identification. To avoid 
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Table 5-1. Reference Star Catalog, Sun Pointing Mode 
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false identification, a bright star is selected that is sufficiently sepa- 
rated from neighboring stars, even if dimmer. The spacecraft is commanded 
to rotate about the sun line, look for this star, and use it to initialize 
the attitude filter. This assures correct identification of the star used 
in the first attitude update. This first update significantly improves 
attitude knowledge and identifications of subsequent stars of 5th order 
magnitude and brighter will pose no problems. 

The required rotation about the sun line is the main reason why it was 
recommended that L-D should align the yaw axis with the sun line and not 
the (canted) solar array normal. Rotation about a non-principal axis that 
is not a spacecraft control axis and to which the star tracker lines of 
sight are not symmetric, appears undesirable. 

The bright and isolated star selected here for initializing the filter 
Is indicated by an arrow in Figure 5-4. It is star number 29 in Table 5-1 
and it has the Yale Number 4023. In the filter simulation, the forward 
looking star tracker encounters this star almost iirmediately and conver- 
gence is rapid as shown in Figure 5-5. An enlarged plot of the same data 
is shown in Figure 5-6 where it is seen that the errors in all axes have 
converged to approximately 5 arc-sec (la) at the end of half a revolution 
(0.5 hr) about the sun line. Gyro biases were calibrated to within 
0.0013 deg/hr (la). Both trackers were read in this simulation every 
30 seconds. . If both trackers are read only every ‘300 seconds, the results 
of Figures 5-7 and 5-8 are obtained where convergence to about 10 arc-sec (la) 
is achieved with gyro bias calibration being about the same, i.e., 

0.0013 deg/hr (la). Both performances satisfy the requirements for suc- 
cessful stellar acquisition and the spacecraft can be considered in inertial 
hold mode. 

With a precision stellar inertial attitude reference successfully 
established, spacecraft ephemeris data is used to obtain an earth pointing 
attitude reference. This may be done as follows. Let Aj^^ denote the direc- 
tion cosine matrix describing the spacecraft's attitude relative to the 


5-25 



WW RRCSEC KlO PITCH HRCSEC ROLL RRCSEC 

I O.S 1 0 200 <00 0 200 400 


ATTITUDE ESTIMATION CONVERGENCE 


.L 

L* 


i.i 11 ^ 


... 

h . 

* 

^ 

_ 

J 






Ffgurs 5-5., Attitude Reference Filter Convergence After Initialization from 
Large Attitude Uncertainties (30 Second Update Interval) 
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Figure 5-6.. Attitude Reference Filter Convergence After Initialization 
(Expanded Scale Figure 5-5) 
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Figure 5-8. Attitude Reference Filter Convergence After Initialization 
(Expanded Scale Figure 5-7) 
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the earth pointing orbital reference frame {x^, y^, z^>. Then is the 
desired earth pointing attitude reference. The stellar inertial attitude 
reference system provides the spacecraft’s attitude relative to ECI given 
by the direction cosine matrix the ephemeris provides the attitude of 
the earth pointing j orbital reference frame relative to ECI by the direc- 
tion cosine matrix A t- Hence, the desired earth oriented attitude refer- 
ence is given by 





The above could have also been done in terms of quaternions and quaternion 
algebra, but frequently direction cosine matrices are required in any 
event and it is advantageous to perform the computations as indicated. 
Figure 5-9 shows a functional block diagram of the on-board attitude refer- 
ence system. 

With Aj^^ established, tne spacecraft is commanded to slew away from 
the sun and assume its earth pointing orientation. After attitude settling, 
the spacecraft is in the normal on-orbit fine pointing mode and ready to 
turn on the payload sensors and take data. Aj^^ represents then a very- 
small-angle rotation matriXjand single axis, i.e. , roll, pitch and yaw, 
control errors can be extracted from the appropriate off-diagonal terms of 
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PROGRAM LfO(TAPE3/ IKPUT^OUTPUT# TAPE !• INPUT# TAPE a«0UTPUT# ' 

C0MM0N/CARR/ATMAST(3) #0IRNAST(3)#HMAST#RDISH# ALPHAS# 

1 ARRW^ARRLKANT/RE ARTH 

COMMON/CORB/GAM«AI#GAMHAD#RL 

C 0 H H □ NV CM LV G A n M A S / G A M M A 1 # G A M M A 2 #~ 6 A M M A L # R T # AL A T 1 
COMMON/CUP /OELTAI# D t LTAL#DELTAT#F1NALT 
Ca«H0N/C0UT/Tl(3)/T2(3);St3» 

DIMENSION A<3)',D{3)#C(3) ' " 

F(ARG)«ARG*3.1A15<>265A/1B0* 

CALL INPUT 

^rewind A“ ■ " ’ 

LAN0»2 

' CALL HEADER(TIME#RL#ftT#GAHMAl#GAHMAS#GAMMAU 

“10‘^coNTiNue — ■ : 

AU)-!, 

A(2) -0. 

A(3)«0. " 

CALL R0T<A^“GAMMA5#2#B> 

CALL >ROT(Bf-GAMMAO#3^ S> 

CALL‘“POT(S/(F (90, )-&AHHAl);i#Bi 

CALL ROT(0,GAHMAL#2#C ) 

CALL ROT(C#F(10O.)#1#S) 

I 1 - X “ 

I2"l 

A(l)-ftT 

CALL"ROt( AV-GAMMAl#3f B) ‘ ^ ' 

CALL CITOO(0#GAHHAL»T1) 

IF(T1(3) ,GT.COS(REARTH> ) 11-0 
CALL”R0T(A/-GAHMA?»3# 6r 
CALL CITOO( B»GAMMAL#T2) 

IF<T2{3),GT.CDS(REAPTH) ) 12-0 

IF(I1.E0,1)~'CALL‘ ARRAir(Ti#S#Rl) “ 

IF(I2.E0.1) call ARRAY(T2#S#R2) 

IF(U.EO.O) Rl— 20. 

iFn 2 .eo.o>'R 2 -- 2 (). 

call eclipse (S,P£AR lH»tSUN# IOAY) 

■ CALL lat(gahmai^gammal#gammai,gammao#alati#alAt#ALOng) 

CALL“GlMBAL{Tl,ftEARlH>GU^G12,G13) 

CALL GIMBAL(T2»RCARTH*G21#G22,G23) 

CALL LANOSEA( ALAT#“AL0NG#LAND) 

CALL~OUTPuT{TlME#Rl#ft2»GAHHAl,GAHRA2#“GAflffALTISCrN#I0AV#^ 

1ALAT#ALONG#G11#G12#G21#G22#LANO) 


• Iw 


00520 CALL UP0ATECTlttt:>0AhMA1^6AMMA2».GAMHAL) 

00530 GO TO 10 

005AO 11 CONTINUE 

00550 END 

00560 . SUBROUTTnE fNPOr 

00570 C0MMQN/CAPR/ATMAST(3) »DIRhAST(3l#HHA$T#R01SH#ALPHAS> 

00560 1. ARRW# ARRL#CANTf PEARTH 

00590 COMKON/CdRB/GANHAI>GA.HMAOVRL *■" " 

00600 COMMON/CML/GAMMAS# GAMMA UGAMHA2#GAMMAL#RT>AL ATI 

00610 CCMM0N/CUP/DELTA1^0ELTAL»DELTAT#FINALT 

00620 DIKENS10N"A( 3) 

00630 DATA ATMAST /3+0 . / t DI RMAST /O * , Oi # -1 . RD ISH/ 3 . / * C ANT /37 . 5 / # 

006A0 lALPHAS/2./>0ELTAl/l5./,DELTAL/220.3U#RL/707A./#RT/A.22AE^f# 

00650 2REARTH/6A.2/,GAMnAl/0Wj.GAHMAi;76. gAMMAS/O* /# GAMMA 1/98.2/ 

00660 3/GAMMAD/37.5/^ALATl/<»l. / 

00670 NAMELIST/IM/6AMMAl>GAMMAL#6AMMASf GAMMAI#PEARTH#RL»RT> ATHAST> 

00680 lDIRMA$T>HMA5T,RDlSH>ALPHA’§>ARftW»AftftCiCAMT»betlAt»f INACT; 

00690 ZDELTAl^DELTAL^GAMMADf ALATl 

00700 F {APG)-ARGt3.1^159265A/lS0* 

’00710 READ (3/ IN)' ^ ' 

00720 on 10 1-1^3 

00730 10 Aa)»OIRHAST(I) 

■007A0 CAiruNlTfAVDIRMAST) 

00750 GAMMA2-GAMMA1-130* 

00760 CANT»F(CANT) 

00770" ‘REAPTH-F <REARtH>' ' 

00780 GAMMAS-F(GA«MAS> 

00790 GAMMM«F(GAMMA1) 

• 00800 GAMMA2-F(GAMMA2 > ' 

00810 GAMMAL*F(GAMMAU 

00820 GAMMAI«F(GAMMAU 

00830 0ELTA1-F<0ELTA1) — 

008AO • 0ELTAL»F(DELTAL> 

00850 ALPHAS«F( ALPHAS) 

"00860 GAMMAD-FCGAHHAO) 

00870 ALATl-F(ALATl) 

00880 RETURN 

00890 END 

00900 SUBROUTINE UPOAT E < T IME# GAHHAl# GAHHA2# GAHMAL » 

00910 COMMaN/CUP/DELTAl#OELTAL»OELTAT#FINALT 

00920 IF(TrHr.6E;FINALTrSTaP — — 

00930 TIME-TIME+DELTAT 



GAMMAl-GAMMAl + OEIUT+ORTAI „ __ 

GAhMAL-GAMMAL ♦ OELT AT+DELTAL 
GAMMA2-GAHHA1 - 130.»3. l*»159Z65AAie0t 

RETURN 

END' 

SUBROUTINE OUTPUT ( TIflE» R1 1 RZ^GAMMAl » GAhNA2» GAMHAU 1 SUNi 

1IOAY#ALAT>ALON6>GU»G12>G21>G22#LAND) 

COhHON/COUT/Tl( 3)#TZ (3)#S(3) ‘ 

50 F0RMAT(lX#lX>F8.3,9(lX»F6.1»>AX*Il#3X^tl»2(l)(»F6.1»»3X# 
F (ARG)-ARG*10O. /3.1*il5Q269A 

AL-F(G'A«MAL)‘ 

A1«F(GAMMA1) 

A2«F (GAMMA2) 

ALAT-F(ALAT) 

ALDNG-F<ALONG) 

GU«F(GU> 

G12-F'(G12) i 

G21-F(G21) 

G22-F<G22) 

WPITe(^#50>flHE# AL,Ai>Rl>GXUG12 

1/A2^R2#G21>G22»10AY,ISUN# ALAT^ALONGiLANO 

RETURN 

ENTRV'HEA'OER " 

51 FORMAT { IX# 5X**T I«E*#3X#*LF/0*#2X*+TDRS1A>2X» 
_^1>C7RANS*>2X#*GIMBAL ANGLES*# IX# TDRS2+# 2X» ♦TRANS*# 

22X#*GIMBAL ANGLES*#2X#*OAY*#lX#*SUN*#ZX#*LAT^#3X#*i.b>4G* 
WRITE(^#51) 

RETURN 

END 

SUBROUTINE ARRAY ( T# S # RE SUIT ) 

COMMON/CARR/ATMAST( 3 )#0IRMAST(3)#HHAST#ftDlSH# ALPHAS# 

X ARRW# ARRL#CANT#REARTH 

DIMENSION V(3)#E1(3)#E2(3I#DT(3>#DM13»#0A(5)#A(3»»T(3)# 
TRANSFORM VECTORS TO ARRAY COORDINATES 

ARRa«-AtAN2<S(1)#S( 3> ) 

CALL COTOA(T# ARRAfDI) 

CALL COTOA(DIRHAST# ARRA»DH> 

CALL'COTOA{ATMAST#A(iRA,OA» 

FIND ANTENNA PATTERN DIAMETER 
TEMP — { (H«AST*0M(3) ♦ 0 A ( 3 ) I /0T< 3 ) t 

IF< TEMP.lt. 0#J* GO T0“'300 

ALPHAO-TEMP - RDISH/TAN(REARTH» 


#'3X#»LAN6*)" 


St3) 


f ' — 



01360 RC-ftOISH ♦ AIPHA04T/N(AIPHAS) . 

01370 C DETERMINE EXTENT OF INTERFERENCE 

01360 CAU BASIS<DT>E1»E2) 

0X390 Ml-0 

Oi^OO M2«0 » 

01610 N«25 

01620 DO 10 I«1>N 

‘01630 Al'-'(27*I-N-‘iVy/N 

01660 DO 10 J-1»N 

01650 AJ- (2.+J-N-1, )/N 

‘01 6 60 DO" 11 K‘-'iV3"~ 

01670 11 V(K)«RC*(AI*EltK> + AJ«E2(K)) 

01680 TEMP-O. 

01690' DO 12' K-lVs' 

01500 12 TEMP-TEMP ♦ V{K)+V(K) 

01510 CHECK-SORT(TEMP) 

01520 IF(CHECK.0TtRCr‘60‘ tO'lO ^ 

01530 hl-Ml + 1 

01560 ALPHA— (HMAST*0M{3) ♦ OAO) ♦ VC3))/OTO) 

01550 DO 13’ K-l»2 

01560 13 A(K)-HMAST*DH(K) ♦ DA(K) ♦ V|K) 4 ALPHA40T(KJ 

01570 TEMP-2. ♦ABS(A(U ) 

0156 0 : 1 F ( TE MP . GTVA ft Rw )' ' CBTCTTO 

.01590 TEMP-A(2) 

01600 1F<TEMP.GT*0. ) GO TO 10 

01610 ‘IF(TEMP,LT*-ARftU‘GD TO lO — 

01620 M2-M2 + 1 

01630 .10 CONTINUE 

01660 RESULT- lOOV* (T.-HZ/Hr) 

01650 RETURN 

01660 300 RESULT-100. 

01670 "RETURN" 

01680 END 

01690 SUBROUTINE COTOA( A.PHO» B) 

01700 C0HMON/CARR/ATHAST(3).0IRHASTt37»RRAnyR0‘ISW?tn»KA57 

01710 1 ARP.W#AftRL#CANTfFEAP7H 

01720 DIMENSION A ( 3 B ( 3 ) #C ( 3 ) 

01730’ ^CAU’R0T{A>Rh0.2.C> 

01760 CALL R0T<C>CANT#1#B) 

01750 RETURN 

01760' END 

01770 SUBROUTINE SASIS(A»B>C) 


__DIMENSI0N_A(3)#B<3I>C(5)>0(3IiEU) 

N*0 

TEMP-0. 
oa 9 1-1^3 

“9“T e Kp=re“HP~+"A''m“wn ^ 

TEST-TEMP/10. 

10 N-N+1 

F (N.'6TV3r'6o“tr'?o ; 

DO 11 I-l#3 

11 0<n«0. • 

DTN)-!. 

CALL CR0SS(A>0»E) 

TEMP-0. 

i)0 i'2 ' “ 

12 TEMP-TEMP ♦ E(n*EO) St 

IF(TEMP.LE.TEST) GO TO 10 

’CAa”UNIT(E»8) " ■■ 

call CR0SS(A.B#D) 

CALL UNIT(D/C) - 

RETURN *■ ‘ - • — - - 

70 DISPLAY+HO BASIS FOUND* 

STOP 

• SUBROUTINE ECL1PSE<S#REARTH#1SUN#IDAYJ 
DIMENSION S(3)iA(3) 

A(l)-0. ■ - ■ 

A<2)«0. 

A(3)-l. 

C A L L"DOT ( SVaV X") 

Y-COS(REARTH) 

rsuN-1 _ 

IF(X.GT.Y)lSUN-0 

IF<X.GT.0.)IDAY-0 

ft ETURN 

ENO 

SUBROUTINE LAT(GAMHA1»GAM«AL>GAMMA1/CAMMA0>ALAT1# At AT a ALONG) 

“^OIMENSION A13).B<3) 

FCARG)-ARG*3,1<»159265A/10O* 

G(ARG)-ARG*10O./3.1A159265A 

— -A(l)-i;; ;; 

A<2)-0. 


52200 

52210" 

52220 

52230 

522^0“ 

52250 

02260 

02270" 

02280 

02290 

02300“ 

02310 

02320 

02330" 

023^0 

02350 

02360" 

02370 

02360 

02390" 

02-«t00 

02^10 

02<t20 

02^30 
OZ^AO 
02A50 ■ 
02A60 
02A70 
02480“ 
02490 
02500 
02510 “ 
02520 
02530 
0254(5' 
02550 
02560 
02570" 
025C0 
02590 
02600" 
02610 


A<3)*0. 

CALL ‘rot (A/-GAHWAL# 2» ft) 

CALL P0T(0#<F(9O.)-GAf1MAn#3>A) 
CALL__ROT(Af F(90. )» 2#ft) 

'CALL ROTTb> (ALAti + GAMrtAD+GAhMAn# 3>Ar 
TLAT«AS1NU(3)) 

A(3)-0* 

'cALr''uNitrA7B) “ 

A(2)-0, 

"A"C3 ) - 0 . : 

CALL OOTtB^A^C) 

ALQNG«ACOS(C) 

IF<8(2) .GT.O.') "A'LDNGi^ALD'NG' 
PARAH-<l*00336417)+*2 
ALAT«ATAN(PARAM*TAN(TLAT) ) 

‘RETURN* 

'END 

SUBROUTINE GIHBAL<T>REARTH#G1^G2>G9) 

‘DIHEnSION T(3)f A(3) ' ■ 

F< ARGJ- ARG*3. 14159269 A A10O* 

G3«0* 

CALL UN r r<T# AT~" 

G1-ATAN2( A(2)/A(D) 

G2-ASIN(“A(3) > 

‘TEttP-REARTH - 90’. 

IF(G2.LT,TEHP) G2-TEMP 
RETURN 

"ENO 

SUBROUTINE C 1 TOD < X> GAttNAL^ T ) 

OINENSION X(3)#Y(3)#Z {3) 
COMMON/COftB/GAMNAl,GAHhAO»H ■■■"" 
F(ARG)»ARG*3 1 14 1592654/ 100. 

CALL ROTtX.-GAMMADf 3.Y) 
"CAU"~ft0T(Y.(F<90.)-GAHMAnn7ir 

CALL R0T(Z.GAHMAL^2#Y) 

CALL RaT(Y>F(10O.)#l»Z) 

‘2(3)-Z<3> ■+ ftL* 

CALL UNIT(Z#V) 

RETURN 

"END 

SUBROUTINE UNIT(XfY) 


OrM6NSION^X(3)>y(3) 

TEHP-Oi ' ‘ "■ 

DO 10 1-1^3 

IP TEMP-TEMP + X(r>+XU) 

f E M ?'• SORT ( fe M P ) “ ^ 

00 20 1-1,3 

lP_Ym-X<I)/TEHP 
RETURN ' " 

END 

, SUBROUTINE DaT(X,V,A) 

DIMENSION X(3')VYt3) “ 

A-0* 

DO 10 1-1,3 

'16‘~A-A' <- “x ( i)* VTn ' 

RETURN 
END 

SUBROuTINTROT <Xt,A,N^XO> 
DIMENSION XI(3>,Xa(3» 

IF(N.E0,2) GO TO 10 

IF < No EO. 3 )” G0'"T0“ 20““' 
XO(l)-XIU) - 

X0(2)-xr(2)*C0S(A) 4 XM3)*SlN(i 
X0< 3 ) — X’l < 2 ) *S I N ( Ar“4 "X r ( 3 ) ♦COS ( 
GO TO 30 

10 XG<i)«Xia)+COS(A) XI(3)*SIN(; 

X0(2)-Xl(2r '■ 

X0<3)»XI(1)*SIN(A) + XI{3)4C0S(/ 
GO TO 30 

•20~xo(i>-xm>»'Cosur»“srf2)»siFt(- 

X0(2)— XI (l)^SIN<A» ♦ Xn2»*CGS' 
xo{3)-xm) • 

“30’ CONTINUE 
RETURN 
END 

SUBPOUTlNe"CROSS(A>6,C) 

OIMENSIDH AC3)#B(3)>C(3) 
C<1>-A<2)*B(B> « AO)*B{2) 

*~C(2)-A(3)*IHl)~Aa)*Bt3) 

C(3 )-Aaj#0(2) - At?)4D(l) 
RETURN 

-^NO 



